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SUMMARY 


Presented herein are the results of a study to investigate the 
theoretical potential of a jet'-flap control system for reducing 
the vertical and horizontal non-cancelling helicopter rotor 
blade root shears. A computer simulation describing the jet- 
flap control rotor system was developed to examine the reduction 
of each harmonic of the transmitted shears as a function of 
various rotor and jet parameters, rotor operating conditions 
and rotor configurations. The computer simulation of the air- 
loads included the influences of nonuniform inflow and blade 
elastic motions. (No hub motions were allowed.) The rotor trim 
and total rotor power (including jet compressor power) were also 
determined. 

It was found that all harmonics of the transmitted horizontal and 
vertical shears could be suppressed simultaneously using a single 
jet control. The jet angle control schedule was found to be 
within practical limits with respect to amplitude. In no case 
did the total jet deflection angle required to suppress shears 
and maintain jet-off trim condition of the rotor exceed 60°. Of 
this maximum jet deflection angle of 60°, no more than 28° was 
required for suppression of the shears. The balance of the jet 
angle deflection was required to maintain the jet-off trim 
conditions of forces and moments. 

The blade response plays an important roll in the mechanism of 
shear suppression. In particular, for the rotor studied, torsion 
was an essential and beneficial element in interacting with the 
jet-flap control to achieve shear suppression. The interaction 
was analogous to that of a trim tab in a conventional airfoil- 
mechanical flap system. The results indicate that the jet-flap 
could be employed to control the torsional response of the blade 
and thereby achieve shear suppression. 

The results obtained concerning the additional power required 
for the jet-flap control rotor system were somewhat obscurred by 
the requirement of the study that the trim (forces and moments) 
of the rotor be maintained to identical values for both the jet- 
on and jet-off operation. In all cases, the power required by 
the jet-flap control to maintain trim appeared to be the dominant 
contributor to the total additional power required (rotor shaft 
power plus compressor power) to maintain trim and suppress all 
transmitted horizontal and vertical shears. 
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1.0 INTRODUCTION 


Helicopter rotor vibration reduction is a desirable and elusive goal 
that has been pursued in many ways and for an appreciable length of 
time. Effort expended has been large because the potential benefits 
are large in terms of vehicle performance, pilot/payload performance, 
and machine life (cost). It was the purpose of the study reported 
on here to continue one of the more promising lines of attack - active 
control of oscillating aerodynamic forces to reduce net shaking shear 
loads transmitted to the fuselage. The particular mechanism for 
exercising this control was chosen to be the jet flap. It is important 
to recognize that the transmitted shear problem belongs to the general 
class of rotor aeroelastic phenomena and that the accurate computation 
of the forcing functions and blade response depends strongly on the 
accurate computation of the induced velocity distribution (nonuniform 
inflow). 

The application of the jet-flap principle to helicopter rotors was 
first discussed over fifteen years ago (References 1 and 2), and 
interest has continued to the present (References 3 to 5, 19, 20, 31 
and 32). Also of interest today are the circulation-controlled rotors 
first investigated both experimentally and theoretically in England by 
the National Gas Turbine Establishment (NGTE) (References 6 and 7). As 
considered by the NGTE, a circulation-controlled rotor consists of 
sections such as circular cylinders and thick ellipses with blowing 
from multiple slots in the vicinity of, but not at, the trailing edge. 
The objective is to obtain lift by controlling the circulation about 
these sections through control of the location of the stagnation 
points. In contrast, the jet flap is applied to sections that are 
more typically conventional airfoil shapes and lead to pressure dis- 
tributions which can differ from that of conventional airfoils primarily 
in that the rear stagnation point is fixed at the jet location rather 
than at the trailing edge. In some current applications, e.g., (Refer- 
ences 8 and 9), the distinction between circulation control and the jet 
flap is not always clear cut. 

The jet flap remains attractive for rotor applications because of its 
potential, not only for integrating rotor propulsion and lift require- 
ments efficiently, but also as a control device for achieving higher 
harmonic control (Reference 2 and 3). Such control can be used to im- 
prove rotor performance, especially at high advance ratios (References 
10 to 12), and can also be applied to reduce vibration. Vibration re- 
duction may be accomplished by eliminating, or reducing, the higher 
harmonic root shears that are transmitted to the rotor hub from the 
blades. The studies of References 13 and 14 have indicated that a 
redistribution of aerodynamic blade loads both spanwise and azimuthally 
can reduce the transmitted shears. This redistribution may be achieved 
via several approaches. One of the approaches is to use the rigid body 
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blade pitch control. The theoretical feasibility of this approach was 
demonstrated in References 13 and 14. Lemnios and Smith (Reference 33) 
analytically evaluated the use of a "controllable twist" rotor to achieve 
the desired loads control. 

McCloud and Evans in Reference 3 noted the possible benefits of reducing 
hub shears by higher harmonic control of the jet flap. Various devices 
have been investigated to achieve the desired control of the jet (Refer- 
ences 9, 31 and 32). The Honeywell system of Reference 31 employs a 
closed loop fluidic system to obtain the jet control while an open loop 
mechanical system is used by Giravions Dorand. Higher harmonic control 
by pulsing the jet supply pressure is being developed by the Navy (See 
Reference 9). These systems represent some of the approaches currently 
being studied to achieve higher harmonic control of the jet. 

Some potential advantages of the jet flap higher harmonic control over 
the "mechanical" control systems are the possibility of a much simpler 
control system at the hub, the higher frequency response which may be 
achievable with the jet control, and the lower control loads possible 
due to the absence of significant inertia loads in the control system. 

Thus because the jet flap offers a promising means of achieving this 
control, the study reported inerein was undertaken to aid in the 
assessment of the potential. 

The objective of this effort was to investigate the theoretical 
potential of a jet-flap control system for reducing the vertical and 
horizontal transmitted helicopter rotor blade root shears. A computer 
simulation describing a jet-flap control rotor system was used to 
examine the reduction of each harmonic of the transmitted shears as a 
function of the jet parameters, the rotor operating conditions, and 
rotor configuration. 

The rotor systems studied were assumed to be shaft powered rotor's, 
i.e. the jet was to be used primarily to control shears and not to 
power the rotor. Conventional collective and cyclic pitch was employed 
to provide rotor control in both the jet-off and jet-on flight conditions. 
In addition to the conventional pitch controls, one jet flap control at 
the first harmonic was provided to allow maintainence of both force and 
moment trim conditions.' That is the forces and moments obtained for 
the rotor in a given flight condition, jet-off were maintained when the 
jet was turned on. The additional jet-flap control was required because 
the rotor trim (primarily rolling and pitching moments) was upset by the 
jet effects. 

Described in the subsequent sections are the rotor systems analyzed, 
the mathematical and computational models employed for the analysis 
and some of the pertinent results and conclusions. 
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2.0 DESCRIPTION OF ROTORS ANALYZED 


The choice of the rotor configurations and the corresponding blade 
properties investigated were made recognizing that no typical jet-flap 
rotor exists. The rotor configurations and blade properties selected 
were based upon the following considerations: 

1. jet-flap rotors currently being considered; 

a) Giravions-Dorand (References 2, 3, 4 and 16) 

b) Lockheed rigid rotor (Reference 17) 

c) The circulation control rotor of NSRDC (References 18 and 19) 

d) The Hughes Tool Company studies (Reference 20) 

2. the primary missions for which jet-flap rotors 
might be used; 

a) heavy lift missions 

b) high speed missions 

3. requirements for reduction of transmitted vibratory shears 
and for simplification of control mechanisms 

4. existing rotor configurations 

The choice for the "basic rotor configuration" was: 

gross weight; 44482 N (10,000 lbs.) 
number of blades; 2 

articulation; flapping hinges for each blade-- 
hinge offset 0.369 m (1.21 ft.) 

blade tip radius; 7.38 m (24.2 ft.) 
blade root radius; 0.921 m (3.02 ft.) 
u blade chord (constant); 0.46 (1.51 ft.) 
solidity; 0.0343 
disc loading; 263 Pa (5.5 psf.) 

Ct/<T= 0.12 

The spanwise distributions of blade properties are basically those of 
the UHIA blade. The geometric properties used for the basic rotor blade 
are presented in Figure 1. The corresponding blade mass-elastic proper- 
ties are presented in Figures 2 through 5. In each of the figures, the 
blade property plotted versus nondimens ionalized radial position has 
been non-dimens ionalized by the respective quantity at the blade root. 
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Presented in Figures 6 through 9 are the mode shapes and frequencies 
calculated for the mass-elastic properties assumed. The inflection in 
the mode shapes near the blade root is due to the hub inertia. The 
rotor operating frequency plot presented in Figure 10 is for no vertical 
hub motion (infinite vertical translational impedance). 

The weight classification of 44482 N (10,000 lbs.) was chosen as repre- 
sentative of a large number of ships (relatively current) made by 
various manufacturers. 

The number of blades was chosen as 2 because this number of blades 
transmits the widest frequency range of vibratory forces to the shaft. 

It is also representative of approximately one-half of the total 
number of ships in this weight classification. 

Flap hinge articulation was selected because it represents approxi- 
mately one third of the rotors in the noted wdight classification. 

The second rotor configuration analyzed was a four (4) bladed rotor in 
the 88964 N (20,000 lb.) class. The rotor blade characteristics were 
maintained the same as those for the "basic configuration" (see Figures 
1 through 10). The solidity for this rotor is 0.0683, Ct/<T = 0.12, 
and the disc loading is 527 Pa (11.0 psf). 

The assumed helicopter configuration is presented in Figure 11. Also 
shown in this figure are the trim forces .and moments and the assumed 
sign convention. The parameters of the assumed configuration are: 

2 2 

Fuselage Flat plate drag area; 8.18 m (26.8 ft .) 

Tail rotor moment arm; 7.93 m (26.0 ft.) 

A conventional control system having collective and cyclic degrees of 
freedom was used to achieve "no blowing" trim of the ship. To facili- 
tate comparisons of the jet-on and jet-off results, it was required 
that the rotor trim conditions (i.e. forces and moments) be the same 
with the jet-on as with the jet^off. Thus an additonal rotor control 
is required when the jet is on^- The additional rotor control is 
provided by the IP jet deflection angle. 

(D Turning the jet on (without deflecting it) significantly alters the 
distribution, over the disc, of the blade section lift, moment, and 
drag and thereby the rotor trim forces and moments. With a conventional 
system, at fixed rotor speed (via the throttle), there are only four 
controls available to control the six rotor forces and moments -- two 
must be left "free", they are the rotor pitch and roll moments. That 
is, the pilot accepts the pitch and roll attitudes required to achieve 
the desired side and propulsive forces. If, however, the moments are 
to be controlled in addition to the forces then it is clear that 
additional controls are required. (See Section 5.2.5) 



This study assumed a "pure jet flap", as opposed to a blown flap or a 
circulation control device, to achieve control^. The jet blowing was 
assumed to be distributed over the outboard half of the rotor blade 
span. The detailed mechanics of the ducting, slot height etc. for the 
jet were not considered in this study. 


In a "pure jet flap", the jet is assumed to emanate from a fixed 
position in the vicinity of the trailing edge of the airfoil thereby 
fixing the rear stagnation point at that location. In a "blown flap", 
the jet passes over a short conventional mechanical flap. In a 
"circulation control" device the jet is blown over a round (or 
elliptical) trailing edge thereby allowing the rear stagnation point 
free to adjust to the local flow requirements. 



3.0 PARAMETER VARIATIONS STUDIED 


3.1 Flight Conditions 

The operating conditions for the basic rotor configuration were: 


Rotor tip speed; 


Rotor tip mach no.: 
Advance ratio: 


RSI - 231 m/s (760.4 ft/sec.) 

XL = 3.4 rad/s (300 rpm = 5.0 HZ) 

= 0.7 (9 altitude of 
762 m (2500 ft.) 

= 0.08, 0.20, 0.30 


Forward Flight 
Velocity : 


Vf = 18.5 m/s (60.8 ft/sec.), 

= 46.4 m/s (152.0 ft/sec.), 
= 69.5 m/s (228.0 ft/sec.). 


The second rotor configuration was investigated only at /J = 0.20 


The ship trim attitude (i.e. shaft angle, o( s ) specified at each 
was determined, based on; 

(i) an estimated fuselage drag which acted through the 
ship center of gravity (c.g.) at all /J and which 
was directly proportional to the freestream dynamic 
pressure, q , 

(ii) a c.g. which was assumed to lie on the rotor shaft 
axis, 

(iii) the assumption that the resultant of the fuselage drag 
and weight was coincident with the rotor shaft axis 


The resultant shaft angles, , as a function of advance ratio, 

>U , are presented in Figure 12. 

With the jet-off, the rotor trim forces were constrained, but the 
rotor pitching and rolling moments (PM and RM respectively) and 
the torque were allowed to be free. 
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The constrained trim forces are in the shaft plane and normal 
to it as shown in Figure 11. The constraint values were chosen 
consistent with requirements on the helicopter weight, drag force, 
rotor torque and shaft angle. 

The rotor torque determined the value required on the side force 
through the constraints: 


T = 

rotor 

Fy rotor = 
where; 

T = 

rotor 

F tail 

Fy rotor “ 
^ tail 


F tail 

x 1 

tail 


(3.1.1) 

F 

tail 



(3.1.2) 

rotor 

torque 



force 

generated 

by the tail rotor 


rotor 

side force 

(in shaft plane) 



distance from shaft axis to tail rotor center. 


When, the jet was turned on, not only were the rotor trim forces and 
shaft angle maintained but in addition the rotor PM and RM were 
constrained to the "jet-off" values. The rotor trim moments were 
maintained by allowing an extra control at IP consisting of the jet 
deflection angle. The amplitude and phase of the IP jet angle (X ) 
required to maintain trim at a given /J and jet blowing coefficient 
are a part of the resultant solution. 


3 . 2 Rotor Blade Tors ional Stiffness 

The blade torsional stiffness assumed for the bulk of the cases 
studied was representative of the torsional stiffness found in 
conventional ( i . e . "no blowing") rotor blades. Since relatively 
(with respect to conventional blades) large torsional moments are 
generated by blowing and because rotor blades employing some form 
of blowing are anticipated to be relatively stiff in torsion, 
additional calculations were made with the blade torsional stiff- 
ness increased by a factor of 10. The 1st and 2nd' tors ional mode 
frequencies were thus increased to 527 rad/s (16. 7P) and 1373 rad/s 
(43. 8P) respectively. The increased stiffness blade was investi- 
gated at y ~ 0.20 with and without jet-on. 

3 . 3 Shears Suppressed 

Given the objective of this study, to determine the effectiveness 
of the jet flap for reducing the vertical and horizontal 
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transmitted blade root shears, the questions to be resolved 
were: 


(i) What value should the shears be reduced to? 

(ii) Should both components (i.e. vertical and horizontal) 

of the transmitted shears be suppressed simultaneously? 

It was decided to suppress all harmonics (i.e. from 2nd to the 
11th)® of the transmitted vertical and horizontal shears to 
zero for all cases studied. Several runs were also made to 
determine the effects on adjacent harmonics of shear when only 
one component of shear at a single harmonic was suppressed. One 
run was made in which only the significant blade root shears were 
suppressed (i.e. up to the 6th harmonic). 

Suppression of the transmitted blade vertical and horizontal 
shears required only one jet flap control mode. This situation 
occurs for two reasons -- first, an ch harmonic horizontal 
blade root shear (in the rotating system) appears in the non- 
rotating system at the ± 1 harmonic orders while the vertical 
blade root shear does not experience this frequency shift; 
second, in a multi-bladed rotor, both the horizontal and 
vertical blade root shears in the non-rotating system cancel 
at all harmonics except 

n = m • NB m = 0, 1,2, ... 

where NB is the number of blades. Thus the vertical and hori- 
zontal blade root shears which will be transmitted to the fixed 
system do not occur at the same frequency in the rotating system. 
Therefore at each frequency only one (either the horizontal or 
the vertical blade root shear) has to be suppressed and hence 
only one control mode is required. Although the simulation had 
a second jet flap control mode available, the simultaneous use 
of both control modes was not required during this study. 

3 . 4 Jet Momentum Coefficient 

The spanwise distribution of the jet velocity was assumed constant 
with radius and azimuth. Thus for a constant ^(with^V) the jet 
velocity will be proportional to the rotor speed and advance ratio, 

V, = k 0 • (R& ) - ( 1 +}J ) (3.4.1) 

® The zeroth and first harmonics of the shears are constrained by the 
trim requirements (See Section 4.3) 
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For this study k e was chosen as 1.2 which is. in the lower end of 
the range to be considered a jet flap rather than boundary layer 
control. 


It was also assumed that the spanwise distribution of the internal 
jet duct pressure was constant, arguing that the internal flow 
losses are just balanced by the centrifugal pumping action. 


The local jet momentum coefficient is defined as 


jet momentum/unit span 
C J freestream momentum/ unit span 

. (a?’) 

pb(R XL) x (r 


(3.4.2) 


(3.4.3) 


where 


<lm 7 

dr = 


P " 


(-P hV^ = jet mass flow rate per unit span, kg/s/m 

2 

freestream density, kg/m 

jetstream density, ( 


b = local blade semichord, m 
K = slot height, m 
V = blade azimuthal position, rad 


"V* = "T/R , nondimensional spanwise distance 


Defining 



spanwise distribution of jet mass flow rate per unit 
span nond imens ional ized to the tip value 



jet momentum coefficient at the blade tip based on 
the rotor speed , flSL 


P b (R n .) 1 
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equation (3.4.3) may be written as 
r = r 

3 Tt • (v»+>L5»«y») 1 


(3.4.^) 


Thus both the jet momentum and its spanwise distribution could be 
varied during this study. 


The values of C 


Jt ® 


selected for this study were: 



0.005, 0.01, 0.02, 0.03 


These values were selected as being in the range of blowing co- 
efficients which might be employed for the purpose of suppressing 
shears on a conventionally driven rotor. 

The jet mass flow rate was assumed constant over the span of the 
jet for all of the cases presented herein. 

The relationships between(c j) , pressure ratio 
( vj ) and nozzle height are based on standard thermodynamic 
equations under the assumption of isentropic expansion. Plots of 
the compressor power required ( ), total jet mass flow rate 

required (mj T ) and pressure ratio vs advance ratio are presented 
in Figure 13. The analyses for these quantities are presented in 
Appendix III. Presented in Figure 14 and 15 are typical plots of 
Cj variation with blade radius and azimuth respectively. 

3 . 5 Jet Control Mode 

As described previously the jet was assumed to be distributed over 
some spanwise extent of the rotor and emmanating from a slot at 
the airfoil trailing edge. The jet angle, 7* , with respect to the 
airfoil chord was assumed to be controllable az.imuthally . The 
spanwise variation of 7« was treated as an independent parameter 
of the problem. 

Two independent jet control modes were provided for. Hence, the 
jet control angle may be described by 

2 « 

r<r,M = £ 2_ (x*. CO sn'f'+Z^snni') (3.5.1) 


I? 



where 



<*•>_ . nondimens ional spanwise mode shape 
for "Cj , j = 1,2 


T* c t * 

W»j , ''wj 


cosine, sine components of the "Cj 
at a given harmonic 


mode 


^ harmonic order with respect to the 
rotor speed 


The jet control mode at IP is used as an additional control to 
allow the maintainence of jet-off trim forces and moments. 

For this study only one jet control mode was required (see 
Section 3.3). Hence the simplest mode shape'Pcj ; was selected, 

•Fcj = i-o 

for all cases studied. 

The values of required to suppress shears are treated 

as unknowns in the computer program. 


3.6 Jet Spanwise Extent and Jet Recovery Factor 

Both the jet spanwise extent and the jet thrust recovery factor, 
Tr , are treated as variables in the program. For all cases 
studied Tr was assumed to be 0.5. 

Two different jet spanwise lengths were analyzed. The bulk of 
the cases were done for a jet spanwise extent of 0.5R. The 
location of. the' jet was over the outer 507. of the blade. Several 
cases were analyzed for a shorter jet extent of approximately 
0.167R located | again, over the outermost blade span (see Figure 
lj jet extent, (1) and (2)). 


3 . 7 Synopsis of Cases Run 

A total of 20 cases involving variations of the above parameters 
are presented herein. An identification number sequence was 
specified to keep track of the various cases run: e.g. 


Case No. 1. 20. 03. 03 

run number within the set 
designates tip jet momentum, Cj T 
designates advance ratio, JUL. 


•-•designates rotor configuration 

(1) 2 bladed, 10,000 lb. rotor 

(2) 4 bladed, 20,000 lb. rotor 
Table I is a synopsis of all cases discussed herein. 





4.0 ROTOR AEROELASTIC SIMULATION 


The Rotor Aeroelastic Simulation (RAS) used in this study is 
described briefly in the following subsections. A more detailed 
presentation of the basic method of approach can be found in 
Reference 13, 14 and 21. Appendix I presents the details of the 
modifications to the aerodynamic portion of the model necessary 
to include the jet. 

In addition to the modifications required to include the jet, an 
extensive effort was expended to include the rotor trim forces and 
moments as dependent variables. The effort followed the lead of 
Chang in Reference 22 and extended it to allow constraint of the 
rotor trim as an integral part of the computational procedure. 

The trim equations are presented in Appendix II. 

This versatile simulation may be employed to solve the 
aeroelastic/trim problem in a direct, inverse or any combination 
there-of manner, i.e. 

1) Given certain prescribed control inputs - determine 
the blade airloads, blade responses, and rotor trim 
forces and moments. 

2) Given certain desired rotor trim conditions and blade 
respsonses - determine the airloads and required 
control inputs. 

3) Any combination of the above (1) and (2). 

It was operational mode (2) which was used during this study. 

4 . 1 Description o f Aerodynamic Simulation 

The aerodynamic model of Reference 21 (essentially identical 
to that of Reference 13) was^employed in the current effort. 
Modifications to the aerodynamic model were made to include the 
effects of the jet on the lift, moment and drag of a rotor blade 
section in forward flight. 

The jet effects were included in the unsteady aerodynamics 
of Reference 21 in a quasi-steady fashion. It is believed that 
this quasi-steady representation of the jet flap aerodynamics 
is a reasonable approach because: 

1) of the preliminary nature of this investigation; 

2) the quasi-steady approximations for the jet flap 
are of the same order as the approximations made 
to define the rotor wake; 
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3) it is believed that the quasi-steady jet flap 
aerodynamics will adequately account for the. major 
effects of the jet flap; 

4) the quasi-steady formulation is adapatable to the 
computational procedures employed; 

5) the theoretical unsteady jet flap aerodynamic 
solutions have been challenged as to their validity 
and this controversy remains unsolved. 


Erickson in Reference 23, has shown that the unsteady jet flap 
aerodynamics are composed of the classical von Karman/Sears 
terms ( L 0| L. j L* J , etc.) of Reference 24, plus similar 

terms describing the effects of the jet, and dependent upon the 
airfoil motions, jet momentum coefficient and shape. 

Because of the above, the results of Reference 21 (and hence 
Reference 13) were followed directly with the appropriate 
modifications to account for a quasi-steady model of the jet. 

The quasi-steady approximation for the jet included the 
effects of airfoil pitching rate as an equivalent camber. The 
steady' state camber effects for a jet flapped airfoil were 
derived by Hough in Reference 25 and provided the basis for 
the quasi-steady approximation. 

No attempt was made to incorporate further approximations, 
for the unsteady jet effects because of the questionable 
validity of the unsteady results which left no means for 
evaluating estimates. 

The jet off aerodynamic model which was used can be described 
as being composed of an adaptation of the classical theory for a 
thin airfoil executing unsteady motions in an inviscid flow. 

The adaptation to a rotor system models the rotor blade as a 
series of spanwise segments whose spanwise aerodynamic properties 
are assumed constant. Each segment was modeled as an equivalent 
chordwise distribution of bound, vorticity whose strength distri- 
bution is determined by the chordwise boundary condition of no 
flow through the airfoil. Since the blade is in an unsteady 
(periodically varying) flow, it is continuously shedding a wake 
which varies both radially and azimuthally due to the radial and 
azimuthal variations of the total bound circulation of the blade 
segments . 1 

This rotor blade wake is modeled as a mesh of concentrated line 
vortices whose strengths and thus induced velocities (non-uniform 
inflow) are always consistent with the azimuthal and spanwise var 
iation of the blade total bound circulation and the blade dynamic 
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response. The periodicity of these variations in the blade bound 
circulation (integer multiples of the rotor speed) when coupled 
with the chordwise boundary condition on the blade airfoil results 
in a set of simultaneous equations in the unknown bound circulation 

Application of the unsteady linearized Bernoulli equation provides 
the relation between the chordwise distribution of bound vorticity 

of the airfoil and the resultant pressure distribution. From 
the unsteady chordwise pressure distribution may be obtained 
the unsteady local lift and moment and hence the spanwise and 
azimuthal variations of these quantities for the rotor blade. 

Added to this model are the quasi-steady effects of a thin high 
velocity jet of air emanating from the airfoil trailing edge 
at an arbitrary angle with respect to the airfoil. The jet 
is modeled as a distribution of vortices and doublets along 
the jet centerline which in the limit of a thin-high speed 
jet becomes simply a single sheet of vorticity on the jet 
centerline. For the steady flow case the effect of this jet 
is shown by Spence (Reference 26) to affect the pressure 
distribution (and hence the airfoil load) in a fashion pro- 
portional to the jet momentum coefficient, airfoil angle and 
jet angle. Hough in Reference 25 included the camber effects 
in a similar fashion. Since a quasi-steady model for the jet 
flap airfoil was assumed, the instantaneous motions of the 
airfoil were treated, insofar as their direct effect on the 
pressure distributions, as though they had always existed. 

The indirect effect of the unsteady jet was approximately 
accounted for via the wake model in that changes in the total 

bound circulation due to the jet were reflected as changes in the 
wake vortex strengths and hence modified the impressed velocity 
distribution on the airfoil. , 

Thus the rotor blade section lift, drag and pitching moment are 
computed including the effects of a wake and nonuniform inflow 
which are always consistent with the coupl'd blade motions and 
the jet effects. 

The details of the equations are presented in Appendix I. 
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4.2 ' Description of Blade Structural Model 

To determine the inertia and elastic forces required for the 
aeroelastic problem, a Lagrangian approach was used. Following 
the development of Reference 27, the kinetic and potential 
energy of an elastic twisted rotating beam was written.. A 
dissipation term proportional to the elastic displacement and 
in phase with the elastic velocity was written to account for 
structural damping forces. The elastic displacements were 
then expanded in terms of the normal modes of the blades, and 
the Lagrangian was applied. Since the forcing function of the 
system, i.e., response independent aerodynamic load, is 
periodic with its fundamental frequency being .R. , the steady- 
state response of the linearized system must be periodic of 
fundamental frequency S!L . Hence, both the motions and the 
generalized aerodynamic forcing functions were expanded in 
Fourier series on SL . Then, by the uniqueness theorem for 
trigonometric series, corresponding coefficients must be equal. 

The result is a set of 2tt (where M is the number of modes or 
degrees of freedom retained) simultaneous algebraic equations 
in the Fourier coefficients (of sine and cosine components) of 
the motions for each of p harmonics. The zeroth harmonic has 
only M simultaneous equations, since there are no corresponding 
sine terms in the series. 

The following assumptions were made in defining the structural 
mode 1 : 

(1) Simple beam theory was applicable. 

(2) Airfoil cross sections were symmetrical with respect 

to the major neutral axis; hence, the center of gravity 
lies in this plane. 

(3) Gravity forces were negligible. 

(4) Flight conditions were steady state. 

(5) Damping (structural) was proportional to the elastic 
displacement and in phase with the elastic velocity. 

(6) Blade elastic deformations were small. 

(7) Pitch axis, elastic axis and center of gravity were non- 
coincident. 

The resultant set of linearized equations of blade motion include 
all the necessary mass, elastic, centrifugal, gyroscopic and aero- 
dynamic force terms required to describe 7 fully coupled degrees 
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of freedom of the rotor blade. The equations have been constituted 
so as to be able to handle blade end conditions of flapping, 
teetering, lead-lag and cantilevered in any combination. The 
aerodynamic forces are described in Section 4.1 and Appendix I. 

For the rotor configuration reported herein, flapping (h^) plus 
6 fully coupled elastic degrees of freedom were provided for. The 
elastic degrees of freedom were: 

1) Two flatwise bending modes, h^ 

2) Two chordwise bending modes, H^» ^ 

3) Two torsional modes,0p©2 

A discussion of the complete set of equations to be solved and 
the method of solution are presented in Sections 4.3 and 4.4. 


4.3 Overall Problem Formulation 

The overall problem formulation and solution procedure used for 
this study are based on those developed by Daughaday and reported 
in Reference 13. The aerodynamic simulation was modified to 
account for the influence of the jet flap as described in Section 
4.1 and Appendix I. The problem formulation was changed to include 
the rotor trim parameters and allow them to be treated as dependent 
variables. 

Daughaday's formulation of the rotor aeroelastic response problem 
represents a powerful and versatile tool. As presently formulated 
it allows the amplitude and phase of any two of eleven dependent 

variables of the rotor to be arbitrarily and independently con- 
strained at each harmonic and the solution to be obtained for the 
remaining nine. Furthermore, it utilizes a sophisticated aero- 
dynamic modeling of the blades and their wake which maintains the 
wake strength distribution and non-uniform inflow distribution 
consistent with the blade dynamic response and unsteady airloading. 

The rotor trim parameters were included within the formulation 
as dependent variables (discussed in Appendix IV) as part of 
this effort. Thus the usefullness of the method has been 
expanded by allowing the rotor trim to be constrained (specified) 
if desired. 



The set of eleven dependent variables are, at each harmonic 
(except the first), the seven blade degrees of freedom (Section 
4.2), two components of blade root shear (horizontal and vertical}^ 
and two control variables. The exception at the first harmonic 
is that the rotor trim forces and moments are included as 
dependent variables in place of the blade root shear variables. 
This is possible when it is recognized that the first harmonic 
blade root shears determine the rotor trim parameters and there- 
fore they cannot be independently specified. Furthermore, since 
these trim parameters are also functions of the other dependent 
variables, their equations can be incorporated into the equation 
set at the first harmonic in a manner similar to the blade root 
shear equations. 

The two control variables at each harmonic (except the zeroth 
and the first) are the deflection angles of the two jet flap 
control modes. At the zeroth and first harmonic the first jet 
flap control mode is replaced by the rigid-body pitch control 
mode. 

Because of the periodicity and the linearization of the structural 
representation as mentioned in Section 4.2 (Structural Model), 
the set of equations at each harmonic are independent. That is 
there is no structural/mechanical interharmonic coupling. However 
the aerodynamic interharmonic/intermodel coupling is included and 
enters through the aerodynamic forces on the right-hand-side of 
the equations (Section 4.4 on solution procedure). 

The overall problem formulation is discussed in more detail 
in Appendix IV. 

This formulation can be used to solve the direct rotor 
aeroelastic problem by constraining the rigid pitch control, 
i.e., given the rotor-operating condition and control 
setting, compute the resulting unsteady airloads, blade 
responses, rotor total loads and moments, and performance. 

It can also be used to determine control settings required 
in addition to the above loads and responses by constraining 
the rotor trim parameters. Or it can be used to perform 
either one of the above solutions with constraints impressed 
at any of the higher harmonics on any of the dependent 
variables . 
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Solution Procedure 


4.4 


The solution procedure employed is an iterative one. The 
rationale for this procedure is briefly described below and the 

logic flow depicted in Figure 16. 


At each harmonic of the rotor speed, let the set of equations 
for the unconstrained variables be written as follows in matrix 
notations, 


£ X* * (4.4.1) 

«a 

(4.4.2) 


where F = (D + C) is the matrix of coefficients relating 

the unconstrained variables and the 
forcing functions 

D is the matrix of mass, spring and 
damping coefficients 

C is the matrix of aerodynamic mass, 

spring and damping coefficients, i.e. co- 
efficients of the response dependent air forces 

X. is the vector of unconstrained variables 

IT is the vector of driving forces (forcing 
" functions). 


If the aerodynamic forcing functions, , (forces independent 
of the blade responses, X« ) and the aerodynamic response dependent 
forces, Cx. , could be explicitly written as indicated above then 
a direct solution would be possible as 


- A., — 

X.* * 


(4.4.3) 


However because the "true" aerodynamic forcing function and the 
blade response aerodynamic forces cannot be accurately expressed 
explicitly, they are computed from an aerodynamic simulation in 
the time domain. But the computed blade section airloads, as 
obtained from the aerodynamic simulation, are the sum of the 
aerodynamic forces representing the "true" forcing function 
(independent of the blade response) and the response dependent 
aerodynamic forces. Thus for a given response vector, Xi , 
the aerodynamic simulation yields the following force vector 


<3F; 


(4.4.4) 
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which is the combination of both the forcing function and the 
response airloads, i.e., C exists only implicitly in the 
aerodynamic simulation. 


Thus Equation (4.4.2) must be written 


D X* = 5 (4.4.5) 

It is the above equation which must be solved. 


If an approximation to the response dependent airloads, CX* , 
is written as EX, then the set of equations can be written as 
follows for solution 


( D + E ) X 
or 

where 

The solution would be 

X 


F X 
F 


6F -*EX 

Co+e') 

F''(<5f+ex) . 


(4.4.6) 

(4.4.7) 


(4.4.8) 


It is obvious from Equation (4.4.6) that if £ were equal to C 
then (4.4.8) would be the solution to Equation (4.4.1), however 
this is not the situation because it is not possible to express 
C explicitly. Equations (4.4.6) and (4.4.8) do however suggest 
a way to obtain the solution by successive approximations 
(iteration). 


For a given (initial) approximation X 9 XJ , a new approximation 
to the 5? can be obtained from (4.4.5) and (4.4.6) as 

( D + E ) * (^.,-CXi) + EXC$ (4.4.9) 

X U , * F‘‘ (6F(-h6)<0 (4.4.10) 

This new approximation is then used in the aerodynamic simulation 
to re-evaluate the force vector which is represented by Equation 
(4.4.4) as, 

GF i+ , * (5 -.-CX; +1 ) (4.4.11) 

If this process of successive approximation is continued (as 
depicted in Figure 16) and it converges, that is 

Uu.-*c) -~ ° 
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then the resulting X will be the solution. This is evident 
from Equation (4.4.9) which would then reduce to Equation (4.4.1) 
because would equal ERl . This is the iterative 

solution procedure used. 

The E-Matrix of coefficients representing the explicit approxima- 
tion to the response aerodynamic loads is based on quasi-steady 
aerodynamics. The D-Matrix represents all the non-aerodynamic 
forces. The F-Matrix (F=(D + E)) is pre-computed for each harmonic 
and inverted for use in calculating the successive approximations 
to the response vector, . 

It should be noted that in this procedure the inflow distribu- 
tion is always consistent with the blade response and airloads. 

The interharmonic aerodynamic coupling enters through the aero- 
dynamic simulation. 

Relative to convergence -- while no quantitive measure has been 
found for ascertaining convergence characteristics of a specific 
solution, it has been found that reducing the sensitivity of the 
right-hand-side of Equation (4.4.8) to the response, "H (i.e., 
improving the E approximation to C ) does improve the convergence. 
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5.0 


RESULTS 


Presented in this section are the synopsized results of the 
investigation. It is felt that the information presented here 
answers the main thrust of the study, namely "What are the 
transmitted shears and what is the corresponding jet control 
schedule and power required to suppress these shears?" This 
information is presented in the form of tables of amplitude 
and phase at each harmonic for all (transmitted and non- 
transmitted) shears and jet control angles. Also presented are 
plots of the azimuthal variations of the total jet control 
angle required. 

The total power required is presented for both jet-on and 
jet-off cases. The jet-on power includes estimates for the 
required jet-flap control compressor power. Jet-off results 
are discussed first, followed by the comparable jet-on 
results. Special jet-on cases are then discussed. 

Additional results in the form of blade responses, bending 
moments, lift load distributions, etc. are introduced only 
as needed to help clarify results or presented as representa- 
tive indicators of what happened in general. 

In all cases in which shear suppression was required, the 
shears were suppressed to zero. 


5.1 GENERAL PRESENTATION OF RESULTS 

5.1.1 Jet-Off Rotor Results 

As described in Section 2.0 the two-bladed rotor was investigated 
at three advance ratios {/ - 0.08, 0.2 and 0.3) with jet-off. 
These cases provided the basic data on the values of the shears 
to be suppressed and the rotor trim conditions to be maintained. 
Presented in Tables 2 and 3 are all of the harmonics of the 
vertical and horizontal (respectively) blade root shears at 
each /J . (As used herein "vertical" and "horizontal" are in 
the shaft-axis system). Also presented are the jet-off shears 
for the four-bladed rotor (Case 2.20.00.02) and the torsionally 
stiff blade (Case 1.20.00.50). 

In all results presented, the amplitude and phase referred to are 
defined by the following general expressions for the tr\ ch 
harmonic variation with y of a quantity, S(\y): 



s„ ( ^ ) = Q« CoS ny +• bn S»A r\v/ 

S - ft* to9lAH*«-Yn) 

k^\_Ci* + b^] v ‘ 

V v, - - tan ' 1 Cb«/a„^ 

The zeroth and first harmonics (OP and IP), of the blade root 
shears contribute to the trim of the craft as described in 
Section 4.3 and Appendix IV. For the two-bladed rotor, above 
the first harmonic,only the even harmonics of vertical blade 
root shears are transmitted, while only the odd harmonics of 
the horizontal blade root shears are transmitted. 

The nominal acceleration level considered to be acceptable by 
pilots at 3P has been found to be about 0.015 g. This value is 
used herein as an approximate criterion of "acceptable" trans- 
mitted acceleration levels at all harmonics for purposes of the 
discussion only. For the investigation, the transmitted shears 
were all suppressed to zero. 

For the two-bladed rotor at /J = 0.08, only the second and fourth 
harmonic transmitted vertical shears exceed the above criterion 
of 0.015 g while none of the horizontal transmitted shears (at 
3p, 5P, 7P, 9P and IIP) exceed this level. At // = 0.20 and 
0.30 only the 2P vertical transmitted shear exceeds the 
criterion and, again, none of the horizontal transmitted 
shears exceed it. 

The effect of increasing the torsional stiffness of the blades 
on the transmitted shears was, generally, to decrease them. 

The exceptions occurred at 4P vertical where the shears increased 
by approximately 57,, and at 7P horizontal where the shear in- 
creased by approximately a factor of 5. In both cases the 
increased values were still below the criterion. 

The rotor configuration change from a two-bladed to a four- 
bladed rotor and increasing the weight from 44482 N to 88964 N 
kept the lift load per blade the same while increasing the disk 
loading and hence the mean downwash. The amplitude and frequency 
content of the downwash was correspondingly altered. Generally, 
the effect on the vertical and horizontal transmitted shears was 
to reduce (or exhibit no change in) their value. The exceptions 
were at 4P and 6P vertical and 7P horizontal where large in- 
creases are noted. 

The changes in the IP shears for the above discussed cases are • 
reflected as changes primarily in the rotor rolling and pitching 
moments. These rotor trim moments are plotted vs u and 


29 



presented in Figure 17. These are the values to which the jet-on 
cases were constrained. Presented in Figure 18 are the corre- 
sponding jet-off rotor power requirements. A minimum in rotor 
power required is noted at about p = 0.125. 

The required control angles (collective and cyclic) are pre- 
sented in Figure 19. Note the slightly larger values of control 
angles required for the torsionally stiffer blade. This is the 
result, at least in part, of the smaller contribution of the 
torsional blade response to the effective angle of attack; thus 
since the rotor is constrained to maintain trim forces in both 
cases, increases are required in the control angles. Note that 
for the case in hand the elastic axis ( e ) is aft of the quarter 
chord (aerodynamic center for jet-off) over most of the blade span 
(See Figure 1) and hence the aerodynamic pitching moment with 
respect to the <2. is nose up . (pos i tive ) . 

A slight decrease in the power required was noted for the 
torsionally stiff blade. 


5.1.2 Jet-On Transmitted Shears Suppressed Results 

5. 1.2.1 Basic Rotor Configuration 

5 . 1 . 2 . 1 . 1 Advance Ratio Variation 

Presented in Tables (4a-c) and (5a-c) are the jet-on vertical 
and horizontal (respectively) blade root shears at all harmonics 
<■ (transmitted and nontransmitted ) for each case run. Presented 
in Table (5a-c) are the corresponding data for the jet angle 
required. See Table I for the case designation. 

The azimuthal variations of the required jet angles for each 
case presented in the Tables are presented in Figures 20 to 26. 

The transmitted vertical shears for the two-bladed rotor which were to be 
suppressed at each /J are presented in Figure 27. (Cases 1.08.00.03, 
1.20.00.09 and 1.30.00.04). The total amplitude (from Table 4) 
at 2P, 4P, 6P, 8P and 10P only is plotted. Similar plots are 
presented for the transmitted horizontal shears (at 3P, 5P, 7P, 9P 
and IIP) in Figure 28. These are the values which were to be 
suppressed. Also presented are the IP shears. which determine the 
trim of the helicopter -- they are termed "TRIM SHEARS". 

The transmitted shears presented in Figures 27 and 28 were 
suppressed to zero at p = 0.08 and p = 0.20 (See Tables 4 and 5). 

The values of the jet control amplitude required at each harmonic 
to achieve this suppression are presented in Table 6 for several 
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values of Cy Ta • The azimuthal variation of the total jet 
control angle (V) at the tip required to suppress these shears 
is presented in Figures 20, 21 and 22. As can be seen from these 
plots, the maximum value of % required to suppress all the 
transmitted shears and maintain the trim forces and moments to 
the jet-off values never exceeded 60°. Furthermore, from Table 
6, the principal harmonic contributors were the IP and 2P 
components. Trim requirements alone always produced large IP 
values as is evident in Figure 22. (See Section 5.2.5) 

It was observed that for the vertical' shears, the trend of the 
required to suppress vertical shears at a given harmonic 
followed the trend of the vertical shear with . Thus, for 

example, the trend of the vertical shears vs Jj for the 4P 

harmonic indicates a decrease in amplitude from /J = 0.08 to 
fj = 0.20, hence the 's required would exhibit the same 
trend. This is true for all harmonics except at 2P where it 

is true only for Cj Tu = 0.01. No such pattern was observed for 

the horizontal transmitted shears. 

No results were obtained with the jet-on at an advance ratio of 
0. 3 because it was not possible to obtain a converged solution 
for this operating condition. The divergence was dominated by 
the 7P response of the second chordwise bending mode which has 
natural frequency of 7.45 P at the rotor operating speed (see 
Figure 10). This problem could be a numerical/computational 
instability of the iterative procedure or it could be indicative 
of proximity to an unstable operating condition for the rotor 
system simulated. It is noted that with the jet-off, at // =0.3, 
a converged solution was obtained but a large region of blade 
stall was encountered on the retreating blade. This could also 
be a significant ingredient in the instability of the solution . 
for this operating condition when the jet is on. 

A method to improve convergence was applied to this case which 
had in similar cases allowed converged results to be obtained. 

The technique involved the determination of correction elements 
to be used in the offending degree of freedom and harmonic which 
aided only the convergence characteristics of the iterative 
solution technique but did not affect the solution once conver- 
gence was obtained. This method was unsuccessful in the present 
cases . 
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1 . 2 . 1.2 


Jet Momentum Coefficient Variation 


Several values of C j t| . were studied at/>-= 0.08 and 
= 0.20 (See Table 1). The azimuthal variation of 
the required "C to suppress shears at juu= 0.08 indicates 
a larger peak ~C for the lower blowing coefficient (Cj r> ). 
However at Jut = 0.20 the opposite is observed. The princi- 
pal harmonic contributions to 'C at all jjl. are the 
IP and 2P components as can be seen in Figure 29. In Figure 
29 are plotted the jet control amplitudes required at each 
harmonic to suppress the transmitted shear at that harmonic 
versus jet tip momentum coefficient. 

No obvious trends with C j T( can be observed except at IP. 
There, for bothjUL = 0.08 and 0.20, the amplitude of the 
jet control angle required increases with increasing blowing 
coefficient values. As will be discussed more fully sub- 
sequently (See Section 5.2.5), this trend is attributable 
to the fact that the lower Cj To , upset trim (i.e. , rolling 
and pitching moments) less from the jet-off values and hence 
require smaller X. to restore trim to the jet off values. 

If only Cj r(> = 0.01 is considered^ Figure 29 shows that at 
all harmonics the required X> trend with jul. follows the 
corresponding flapwise or chordwise shears trend with jLL. 
as pointed out in Section (5. 1.2. 1.1). However for increased 
C = 0.03 it is observed that these trends are reversed at 
2P, 3P, 5P and 7P. These results provide some small indica- 
tion as to the complex nature of the interaction which occur 
between the influence of the jet and the blade responses. 

The effects of blade resonances, interharmonic couplings and 
the manner in which shear suppression is achieved via the 
balancing of inertia, spring, damping and aerodynamic forces 
change with C-y T and ju- . Though not definitive, a dis- 
cussion of these interelations is presented in Sec t ion 5 . 2 . 3 . 

Presented in Figure 30, as indicative of what happens to blade 
responses when the jet is activated, is the azimuthal varia- 
tion of the tip deflection of the 6 elastic blade degrees 
of freedom for the jet off and Cj Tt = 0.03 case a,t 
M- = 0.20. Note in particular the greatly increased 
torsional response of the blade with the jet on. Both 
first and second torsion are affected (note scale change 
between 1st and 2nd torsion). If the variation of first 
torsion in Figure 30 is compared to the T. of Figure 22 . 
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for all shears suppressed (the torsion and t both 
result from Case 1.20.03.08) the following are noted: 

1) a very strong similarity in the azimuthal 
variation between 8, and "C 

2) approximately a factor of 10 difference 
in the ©, and X peak amplitudes 

3) an approximate 180° shift in phase 
between 6, and “C .. 

5. 1.2. 1.3 Torsionally Stiff Blade 

The observations in the preceeding section concerning 
the blade torsional response led to the speculation 
that a torsionally stiffer blade would require larger 
jet deflection angles, X , to accomplish the shear 
suppression. The results for a torsionally stiff blade 
(Case 1.20.03.58) are presented in Tables 4, 5 and 6 
and in Figure 23. The torsionally stiff blade did 
require significantly more,'t , to suppress the shears 
at all harmonics except 2P and 8P (see Table 6). 
However, the stiff blade required' less X at IP for 
trim but the IP is still the dominant frequency com- 
ponent. A discussion of these results is presented 
in Sections 5.2.1 and 5.2.2. 


5. 1.2. 1.4 Selected Shear Suppression 

Several cases were run to determine the effects on the 
transmitted shears when only selected shears were suppressed. 
The cases considered were at /-> = 0.20 for Cxf 0 ~ 0.03. The 
shears chosen for suppression were the 2P vertical and the 
5P vertical shears. Also run was a case in which the 2P, 

4P and 6P vertical shears and the 3P and 5p horizontal shears 
were suppressed (i.e., all transmitted shears through the 6P) 
In all cases, trim was maintained at the jet-off values. The 
required X azimuthal variation for these cases is presented 
in Figure 25. Also presented in that figure is the X for 
suppressing all transmitted shears through the IIP harmonic. 
The X variation required for suppressing all transmitted 
shears compared with that required for suppressing only the 
first six transmitted shears can be seen to differ mostly at 
the 7P and 9P components. These 7P and 9P components of X 
are those required to suppress the horizontal shears at 
7P and 9P. 
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Several interesting observations can be made from these 
cases, namely; 

1) interharmonic effects may be quite strong 

2) suppression of one harmonic of shear 
(flapwise or chordwise) will affect the 
magnitude of all others -- generally the 
shears at harmonics immediately above and 
below the suppressed harmonics are most 
affected . 

A quick scan down the vertical and horizontal shears of 
Tables 2, 4 and 5 for Cases 1.20.00.09, 1.20.03.09 and 
1.20.03.10 will verify the above. Also observe Case 
1.20.03.20 where only rotor trim was constrained; note 
especially the very large change in the 2P components 
and the substantial changes in the chordwise shears 
induced by the jet at the 6P and 7P harmonics. Thus 
simply turning on the jet and using only a IP variation 
to maintain trim caused large changes in blade root 
shears at other harmonics. It is speculated that the 
large changes at the 2P and 7P harmonic are at least in 
part related to the fact that the chordwise bending 

mode resonances are in these ranges (See Figure 10 ). 


5 . 1 . 2 . 1 . 5 The Short Span Jet 

To determine the effect of changing the spanwise extent 
of the jet a "short jet" case (1;20.03.S6) was run 
suppressing only the 5P harmonic of the vertical shear. 

This harmonic was chosen because it was closest to the 
first torsional natural frequency. When compared to 
the long jet (See Figures 25 and 26 and Table 6) -- 5P 
shear suppression (1.20.03.10), a smaller IP variation 
to maintain trim was required by the short jet; however 
almost twice the T amplitude at 5P was required to 
suppress the vertical shear. The more important aspect 
of this case is the fact that the power is much reduced 
for the short jet (See Section 5.1.3) 

5. 1.2. 2 Four Bladed Rotor Configuration 

The four-bladed, 88964 N rotor results are presented in 
Figure 24 and Tables 4 through 6 to Case 2.20.03.02. All 
transmitted shears were reduced to zero. Relative to the 
two-blade rotor, large changes in X requirements occurred at 
virtually all harmonics; except for 2P all even harmonic 
requirements were substantially increased. No pattern for X 
required to suppress the horizontal shears emerged. 
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5.1.3. Total Power Considerations 


As discussed previously, the analysis computes the basic 
rotor power required jet»on or jet-off due to the aerodynamics 
and inertial loads of the blades. To this power was added 
the estimated compressor power required to supply the jet 
(See Section 3.4 and Append ix III for details). Thus the 
total power required was given by 


Tt * P R +P^ 


(5.1.3. 1) 


The compressor power P c was derived as a function of 

and jul and is presented in Figure 13. The rotor shaft power, P R , 

was determined as 


P B = 



2n 

* 

R 

** 

< h 



<£* <*l +£) 







(5.1.3. 2).. 


'**- 
2 = 
2? = 
"i* = 
M = 
SL = 


power supplied to the shaft (SHP) 

blade section lift including jet contributions 

blade section drag including jet contributions 

blade section induced angle 

number of blades 

rotor speed 


(See Appendix I, Equations 1-31 and 1-37 and Appendix III, 
Equation III-2) Recognizing that equation (5.1.3. 2) can 
be written as 


^R a + ~"Pcj 


(5.1.3. 3) 
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where 



“power due to the profile drag 

=power due to induced effects 
“power (propulsion) due to the jet 


Equation (5.1.3. 1) can be then written as 


■MlVP* -h T ) 


Recognizing that 

oc 

Tx oc j£-oll 

CC 


(5.1.3. 4) 


and that , o(; , T and may be represented by a 
Fourier Series on , it can be shown that the "mean" 
total power, "Pf , i- n the Fourier sense > depends not only 
on the steady valves of oL^ , , "C and o£ but also on 

the higher harmonic amplitudes of these quantities. 

For example, 
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(5.1.3. 5) 
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(5.1.3. 6) 


where 

C | “function of only the steady valves of 0 
“constant of integration 

©^ 9 =cosine, sine components of the n ch harmonic 
of the sum (o^-FC ) 



Thus while ordinarily the jet-off harmonic variation 

of ,aL» and X may be ignored with regard to power, 
with the jet-on they may become important. Further it 
should be noted that there is a power penalty to be paid 
for activating X which is proportional only to the jet 
deflection amplitude required regardless of the harmonic. 

Also keep in mind that C$, o and X also influence X an ^ 
oLi . Hence the overall interaction is quite complex as 
discussed in the following paragraphs. 

The total power, P T , required at y = 0.20 to maintain trim 
and suppress all transmitted shears is presented in Figure 
31 as a function of C Jr/ • It is observed that this total 
power is a monotonical ly increasing function of Cj To , 
i.e., in these results, there is no indication of a minimum 
with respect to . However these results do not 

necessarily rule out the possibility of such a minimum. 

The lack of a minimum with respect to ^To may be the 
result of the trim philosophy adopted for this study (see 
Section 5.2.5). 

The following observations can be made concerning the 
various contributions to the total power required. It is 
helpful if Equation (5.3.1 4) is rewritten as the following 
two groups of terms, 

P T = p.(l- P<_j/P t ) r + Pi) (5.1.3 7) 

Concerning the first group of terms, it is first recalled 
that the compressor power required, P t , is a monotonically 
increasing function of Cj T j> as in Figure 13, and that the 
jet reaction/recovery power, P Cj , is that portion of P L 
recovered (i.e., Pc* / Pc < 1 ). Second, it is recalled 

(Equation 5.3.1 5 and 6) that P Cj increases with Cj- g 

and decreases with (regardless of the harmonic) 

the required t's , however, depend on Cj Ti# . That is the 
Z 's required to suppress the shears, , decrease with 
increasing Cy Tg and the X 's required to maintain the 
jet-off trim, Z T , increase with C *t 8 . Thus , if only the 
X required for trim, Xr , is considered, then the net 
power resulting from the first group of terms is a mono- 
tonically increasing function o It is not possible, 

apriori, to determine the variation of this net power with 
C.t To due to the X 's required for shear suppression, X$ , 
because it will depend on the relative rate of change of 2T S 
and ’Pj. with 
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However, for the trim philosophy used, the increase of 
with Cjt. is much greater than the decrease of the *ts , 
and furthermore, because p tJ varies as the square of 'C 
(Equation 5.1.3. 6), it is believed that the "C T inf luence 
dominates and thus the total power p T increases mono- 
tonically with Cj Tft 

The influence on the total power represented by the second 
group of terms (Equation 5.1.3. 7) is more obscure than 
the first group of terms. Both the "profile power", P c# , 
and the "induced power", P x , will be indirectly influenced 
by Qt-t and £ . It is expected that C~ Ta would reduce 

* O 

the profile power, P Cfc , but the magnitude of the influence 
would depend on the operating state of the local airfoil 
section (i.e. Reynolds number, Mach number, angle of attack, 
etc.). Th'e influence of Cj To on the induced power, Pj , 
would be through the lift distributions and thus the induced 
velocities. However it is believed that it would, to a 
larger extent, be determined by the specific application of 
the jet-flap. That is, it would depend on the redistribution 
of the lift required (and thus the ~C required) for the 
specific application. 

Also plotted on Figure 31 is the total power, P T , required 
to turn the jet on and only maintain trim at JJ =0.2 and 
Cy To = 0.03; it is, 

P T = 1193 HP. 

Of this total, 872 HP was the compressor power required. 

The jet-off shaft power, P T , was 892 HP. 

Thus a 347o power increase was required just to maintain 
rotor trim at its jet-off value when the jet was turned on. 
Further it is noted that for this case, the compressor, 
power required just about equals the jet-off power required. 
It has been speculated that,if trim had not been required 
to be maintained in the manner of this study, the power 
penalty indicated herein may have been substantially reduced. 

Even with the adverse power requirement (relative to trim 
maintainance ) , Figure 31 implies that the job of suppressing 
all transmitted shears can be accomplished simply by re- 
ducing the jet tip momentum sufficiently. The increment in 
power required (See Figure 31) to go from maintaining trim 
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only to trim plus full shear suppression is indicative 
of the effect of higher harmonic X contributions. An 
additional 117 HP was required just to suppress all 
transmitted shears while maintaining trim (compressor 
power held constant for Cj 7o = 0.03 and /s = 0.2). Thus 
the total power required was increased by 477. over that 
required jet-off, but only 137. of the increased power 
required was attributable to the requirement to suppress 
all transmitted shears. 

What was not definitively ascertained in this study was 
the corresponding limits on the resulting jet control 
angle amplitudes. For the lower, Cj To (0.005) studied 
the total value of X , required to maintain trim and 
suppress all shears was not large (see Figure 21). 

This was due to the fact that as Cj T(> was decreased the 
rotor trim was upset less and less and hence the 1P-T 
required to maintain trim was smaller. Thus because the 
IP— "TT component is by far and away the largest term in 
the series representation of X , the increased values 
of V required at the higher harmonics to suppress shears 
at this low Cj To , were more than compensated for by the 
reduced 1P-T; requirements. 

At some sufficiently small value of Crr* it is expected 
that, even though the 1P-"C required for trim would be 
very small, the required X to suppress higher harmonic 
shears would begin to become prohibitively large from a 
practical standpoint. Once the jet angles become 
sufficiently large, the analysis employed herein begins 
to suffer from certain small angle assumptions made in 
the aerodynamic modeling. 

The effect on power required for suppressing selected 
shears is presented in Table 7, There it is seen that 
for a f ixed // and Cj To ( 0.3 and 0.03 respectively) requiring 
successively more shears to be suppressed increases the 
P T required. The cases suppressing the 2P and 5P flap- 
wise shears tend to substantiate that the power required 
to suppress a given harmonic of shear at a given Cj T 
and p goes like the X required. 

The power required for the torsionally stiffer blade is 
considerably reduced over the power required for the 
less stiff blade (Figure 31). This power difference is 
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attributable to the large changes in the't required and 
the large changes in torsional response. The jet-off 
power requirements for the torsionally stiff blade were 
also found to be slightly lower, 874 HP, than for the 
less stiff blade. Thus stiffening the blade had a bene- 
ficial effect on the power required. 

The reduced IP torsional response is reflected directly 
in the reduced l requirement at IP (see Table 6b and 6c - 
Cases 1.20.03.08 and 1.20.03.58). A large reduction in 
the 2P component is also noted. Although not explicitly 
identifiable, it is believed that the large power reduc- 
tion observed is due to the reduced IP- T. and the reduced 
shear load to be suppressed at 2P because of the reduced 
torsional response. The X requirements at all the remain- 
ing higher harmonics are significatly increased. The 
reduced IP and 2P X. requirements more than outweigh the 
increased higher harmonic X requirements. Thus the over- 
all result is the large reduction in power required. 

These explanations are somewhat speculative because the 
requirement for maintaining trim tends to obscure the 
influence of the other participants. 

Similarly, the comparison of P T between the long jet and 
short jet cases (1.20.03.10 and 1.20.03.56) in Table 7 
indicates a much reduced power requirement. Thus, although 
a much larger X is required (but not excessively large), 
the same trim balance and shear suppression was achieved 
at a much smaller power. 

Presented in Figure 32 is a crossplot of P T vs p for 
the cases run. Power jet-on trends follow the jet-off 
trends as might be expected. 

the power required jet-on and off for the four bladed 
rotor is presented in Table 7 (Cases 2.20.00.02 and 
2.20.03.02). A comparison of the jet-on power for the 
two bladed rotor and the four bladed rotor indicates that 
approximately the same power penalty is paid to suppress 
all shears for both configurations. 
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5.2 Discussion of Results 


The results as presented in Section 5.1 encompass the complete 
scope of this investigation to assess the theoretical potential 
of a jet-flap control system for reducing the vertical and 
horizontal transmitted blade root shears. In this section, 
some of the more pertinent or interesting aspects of these 
results are noted and discussed. It should be recalled that 
the rotor aeroelastic response problem is not a simple one. 

The many parameters, forces, and variables are intimately inter- 
related, therefore, simple explanations may not always exist. 
Furthermore those which are presented may not necessarily be 
generalized beyond the specific problem in-hand. 


5.2.1 The Role of the Blade Torsional Response 

V. 

The results of some of the inital calculations for suppression 
of the transmitted blade root shears on the basic rotor con- 
figuration, suggested that the blade torsional response 
may be a primary element in the jet flap control system. It 
was observed that because of the effectiveness of in 
generating blade section aerodynamic pitching moments (relative 
to the elastic axis), it was the dominant driving "force" on the 
torsion response. That is, tl was controlling the first torsion 
mode. 

The question was then raised as to which was the "effective" 
aerodynamic force control in reducing the shears--the jet 
deflection angle, “t , or the blade torsional deflection angle, 

0 ? Actually TJ is the "primary" control (because it is the 
input to the rotor) and both X. and © influence the aero- 
dynamic loads but what are their "effective" aerodynamic roles? 

Is the torsional response reducing the effectiveness of 'C 
or is it a beneficial influence on the efficiency of the jet 
flap control system relative to shear reduction? With the 
objective of obtaining information relative to the above kinds 
of questions, additional calculations were made. 

In order to assess the importance of the torsional response 
. relative to X , the torsion response was virtually eliminated 
by increasing the generalized stiffness in torsion by a factor 
of 10 (so that @ 16. 7P and (9 43. 8P). The X control 

schedule required to suppress all the transmitted blade root 
shears was then re-calculated. The resulting azimuthal variation 
of V is presented in Figure 23 and the corresponding harmonics in 
Table 6c (Case 1.20.03.58). The harmonics of the resulting blade 
root shears are presented in Table 5c. The corresponding results 
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for the standard blade are in Figure 22 and Tables 6b and 5b, 
(Case 1.20.03.08). It is observed that, virtually eliminating 
the blade torsional response has increased the X -amplitud e 
requirement at all harmonics except 1, 2 and 8. This would 
seem to indicate that the torsional response of the standard 
blade is generally beneficial relative to the effectiveness of 
the jet flap in shear suppression. However it degrades the 
effectiveness relative to trim at IP as described below. This 
result further supports the conjecture that (for the standard 
blade) the jet flap is controlling torsion and the torsion is 
the "effective" aerodynamic force control. That is the jet flap 
is operating as a control tab (trim tab) similar to the Kaman 
"controllable-twist" rotor (Reference 33 ). 

The above exception at the 1st and 2nd harmonic (¥\ = 1 and 2) 
are not really understood. However it should be noted that the 
constraints are different atn=l, i.e., the rotor trim force 
and moment are constrained rather than the blade root shears. 

It may be that, relative to rotor trim,<ti s the dominant influ- 
ence and the torsion influence is opposing that of X , i.e. 
reducing its effectiveness. It is noted that, at R = 1, 
torsion is nearly 180 degrees out of phase with X (i.e., © 
is negative when X ls positive). Thus, for the stiff 
blade at fl = 1, the reduced torsion response (from 4.8° to 
0.4°) would increase the effectiveness of X relative 
to trim and thereby explain the reduced X requirement (from 
31.1° to 27.4°). No rationalization of the anomalous results 
relative to the X requirement at h = 2 and W = 8 can be 
offered. 


The greater effectiveness of the blade section angle-of-attack 
© , over that of the jet-deflect ion angle, X , in controlling 
the aerodynamic forces arises because the X- dependent force 
coefficient slopes (e.g., Cu^. ) go to zero asCt^ goes to zero 
while the ©-dependent slopes (e.g., Gi* ) go to their jet-off 
values (see Appendix I). Thus for example, the variation of 
the ratio, C ttl / » (the respective at and X lift curve 
slopes) with gj is as sketched below. In the low range of 
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blowing coefficient, of interest herein, it is observed that the 
effectiveness of 6 is very much greater than X . The situation 
is similar for the blade section .drag and moment. 

The above results suggest the possibility of a' l jet-flap torsionally 
controlled rotor” (JF/TCR) , similar to the Kaman CTR, that is, a 
rotor designed so the primary function of the jet flap is to "drive" 
the blade torsional response at each harmonic. The aerodynamic 
force control would be obtained from the angle of attack variation 
due to the torsional response at each harmonic. The efficiency 
of the jet-flap control system may be improved by this mode of 
operation. For a given blowing coefficient reducing the span 
of the jet will reduce the total mass flow rate required and 
thus the compressor power. By locating the jet at the blade tip 
where the torsion deflection is maximum, its effectiveness 
(generalized force in torsion mode) is maximized. 



5.2.2 Jet Flap Torsionally Controlled Rotor 

To obtain an initial assessment of the potential of the 
JF/TCR described above, the span of the jet was 
reduced by two-thirds. This configuration was then used 
to repeat a previous case (1.20.03.10) where only the 
dominant 5th harmonic vertical blade root shear is suppressed. 

The required azimuthal variation of 'C is presented in 
Figure 26 and the corresponding harmonics in Table 6C 
(Case 1.20.03.S6). The corresponding blade root shears 
are presented in Tables 4C. and 5C. The comparative results 
for the basic configuration are presented in Figure 25, and 
Tables 6B, 4B and 5B (Case 1.20.03.10). 

The results are very encouraging for this case, the power 

required was reduced from 1200 HP to 892 HP and, while the 
required 't at the fifth harmonic increased from 3.4 degrees 
to 7.3 degrees, the torsion response only increased 20 per- 
cent (from 1.23 degrees to 1.48 degrees). This result further 
illustrates that the torsion response is the fundamental aero- 
dynamic control and . ‘tl is controlling torsion as a control 
flap. The first harmonic component of required to maintain 
the jet-off trim was reduced from 28.2 degrees to 24.5 degrees. 
Apparently the shorter span jet upsets the trim less and thus 
requires less first harmonic control to maintain the jet-off 
trim (see Section 5.2.5). 


5.2.3 Mechanism of Shear Suppression 

At each harmonic of the rotor speed, there are spanwise 
distributions of aerodynamic and inertia loads proportional to 
the blade response in each of its degrees of freedom. There is 
also a spanwise distribution of the aerodynamic loads which are 
independent of the blade response (i.e., the aerodynamic 
forcing function). The blade root shears (vertical and 
inplane) are the sum of these spanwise distributions of 
aerodynamic and inertia loadings. 

There are conceptually two possible ways of eliminating a 
given harmonic component of shear. The first is, to 
effectively cancel the aerodynamic forcing function in each 
blade degree of freedom — this would obviously reduce all 
the blade response at that harmonic to zero and thereby 
eliminate the response dependent aerodynamic and inertia 
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loadings in addition to the forcing function airloading. The 
number of independent control modes required for this con- 
ceptual approach would have to equal the number of blade-degrees - 
of-freedom -- this seems impractical at the present. 


The second obvious and more practical way of eliminating a 
given harmonic component shear is to only require that the 
aerodynamic and inertial contributions sum to zero. This 
requires only one independent control for each component 
of shear to be constrained at each harmonic and is the method 
used in this study. 


Only as an example to illustrate how the various contributions 
to a blade root shear can change to make sum zero at a given 
harmonic, these components are presented in the vector plots, of 
Figure 33 for case (1.20. 03. S6) of suppressing the dominant 
5th harmonic vertical blade root shear with the shortened jet span 
Individual aerodynamic contributions for each blade mode and 
the forcing function are not available separately (as explained 
in Appendix IV) but only their sum. Thus presented in Figure 33A 
are all the significant contributions to the 5th harmonic blade 
root shear before it was suppressed, i.e., with the jet off 
(Case 1.20.00.09) -- they are: 


A - 
V 

3V 


the total aerodynamic component 
flapping inertia load 

1st flapwise bending mode inertia load 
2nd flapwise bending mode inertia- load 
1st torsion mode inertia load. 


When this 5th harmonic shear is suppressed with the short jet 
the components re-adjust as shown in Figure 33 b. With the 
jet off, the largest contribution to this shear was the 2nd 
flapwise bending mode inertial load -- this mode is near 
resonance (flJ»» = 5.27). When this shear is suppressed to zero, 
the components change as follows: 
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Thus, all the components increased except that due to the 
bending response near resonance. 

This general increase in the magnitude of the individual 
components of the reduced total shear implies that the re- 
sulting total shear could be sensitive to the amplitude and 
phase of the jet control angles. That is, open-loop opera- 
tion may be unsatisfactory because it could result in an 
increase in the level of the shears if the jet control angles 
were not properly chosen. 


5 . 2 . A Blade Dynamic Response to Shear Suppression 

. = As was illustrated in the previous section, when the blade 
root shears are constrained, the blade dynamic response in 
all modes is required to re-ad iust so that the appropriate 
aerodynamic and inertia contributions sum up to the constrained 
values. In general , this resulted in an increase in the blade 
bending responses at frequencies away from resonance and a 
reduction of these responses at frequencies nearest resonance. 

This can be observed in Figure 30 where the azimuthal varia- 
tions (time histories) of the response in each blade mode are 
presented for the jet off case, and for the case with all the 
transmitted shears suppressed. It is observed that the 1st 
flap-bending mode 2.94) response at A. =3 is reduced, the 

2nd flap-bending (6J n = 5.38) response at A = 5 is about the 
same, the 1st edgewise-bending (UJn=0.98) response atn =1 is re- 
duced, but the 2nd edgewise-bending (W*\ = 7.45) response increased. 
As discussed in Section 5.2.1, the blade torsion is strongly 
coupled to the jet deflection angle, X » so that the jet 
control drives the torsion response. Therefore, as observed 
in Figure 30 the blade response in both torsion modes is 
increased when the jet flap control is used to suppress all 
the transmitted shears. 



5.2.5 Jet Influence on the Rotor Tr im 

In this study the tip jet momentum coefficient, is 

constant with azimuth, vp , however, as explained in Appendix. 

I, the blade section jet momentum coefficient, Cj , varies 
through relatively large amplitudes as depicted in Figure 15. 
Because of this large variation of Gj with sy there is a 
significant redistribution of the blade airloads over the 
. disk and a resulting change in the rotor total forces and 

moments (especially the moments) when the jet is turned on, 
i.e. the rotor trim is upset. 

With' the jet off, the three rotor trim forces are controlled 
by the three components of blade rigid body pitch, i.e., the 
collective pitch, lateral, and longitudinal cyclic pitch. 

The yawing moment is controlled by the "pedals" (tail rotor 
collective). The resulting rotor pitch and roll moments are 
accepted and thus determine the rotor shaft orientation 
(attitude). . With the jet on, the philosophy adopted for this_ 
investigation was that the rotor ‘forces and moments shall be 
the same as for the jet off. 

To control the rotor moments in addition to the rotor forces 
requires two additional controls -- in this study, they are 
the 1st harmonic sine and cosine components of the' jet 
deflection angle, t ( see Section 4.3 and Appendix IV). 

This requirement, for the rotor moments to remain unchanged, 
results in the same rotor shaft angles (i.e., fuselage attitude) 
with the jet on as with the jet off. This trim procedure is not 
necessary -- the resulting, different, jet-on moments could 
alternatively have been accepted as an attitude change. 

It should be noted that the first harmonic blade pitch and 
jet angle control requirements are uniquely determined for a 
specific rotor trim condition. That is they are not inter- 
changeable when both the forces and the moments of the rotor 
are specified because the jet angle control results in a 
different distribution of blade lift, drag, and moment over 
the disk than does the blade pitch control. However if only 
the rotor forces (or only the moments) were specified then 
either control could be used independently or they could be 
used on arbitrary combinations. Thus the specified rotor 
forces (or moments) can be obtained with any combination of 
pitch and jet control but then unconstrained forces (or 
moments) would be different for each combination. 
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Because the amplitude of the variation of Cj , 
with V (Figure 15) increases with increasing , the 

change in the rotor trim due to the jet increases and there- 
fore the amplitude of 1st harmonic V required to restore 
trim also increases with C 3r . This variation (with Cj lo ) 
of the 1st harmonic X. required for trim is presented in 
Figure 22. It is observed that amplitude of X required 
at l\ = 1 to maintain the jet off trim is quite large, varying 
from 13.6 degrees at C Jr. = 0.005 to 31.1 degrees at 
Cj r = 0.03. These relatively large jet deflections for 
trim result in signifcant power penalties as explained below. 
They are also believed to be the source of the relatively 
large X requirement generally required at Y\ =2 (via the 
inter- harmonic aerodynamic coupling). 

When the thrust recovery factor, IJi > is less than one, the 
jet provides a steady propulsive torque to the rotor which is 
a function of the total angle of the jet relative to the shaft 
plane. (A value o f T r = 0.5 was used for this study). This 
total angle is the sum of the blade section geometric angle 
relative to the shaft plane and the jet deflection angle,!; , 
relative to the blade section chord. As a result of the inter- 
harmonic coupling, when the total jet angle, cK. n , at each har- 
monic is less than one radian, the steady propulsive torque, 
Qr, derived from the jet decreases approximately as the 
square of each harmonic of the jet angle, 0„ , i.e., 

Gv«{l - + I v, 9J .e" } 

That is, the rate of decrease of Q* increases with the larger 
angles. Thus the relatively large X required to maintain the 
jet off trim resulted in large power increments -- for example, 
at y =0.2 and C;j T# = 0.03 this power increment was about 
200 horsepower. This is the total increment in power 
which includes the profile and induced components as well 
as the Q^. component. 

These results, relative to the "trim-power", raise the follow- 
ing question. Would the total power required with the jet on 
be less if the rotor total moments were accepted an an attitude 
change? A proper answer to this question, of course, would be 
dependent on the drag variation, with attitude, of the specific 
fuselage considered. This is the main reason, for the choice 
made herein, to require the rotor moments to be the same with 
jet on and jet off. However, this alternate trim procedure 
should be investigated, at least for a typical fuselage 
drag-characteristic because the procedure used herein may 
represent an unnecessary power penalty. 
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6.0 CONCLUSIONS 


The general overall conclusion of this study is that the rotor 
jet-flap control system appeard to be( theoretically) a prac- 
tical means of achieving efficient higher harmonic control 
which could be used for many applications. 

The results indicate that a jet flap control system has the 
potential of reducing all of the vertical and inplane trans- 
mitted (non-cancelling) blade root shears simultaneously with 
a single jet control mode. Furthermore the results indicate 
that the control angle schedule and additional power required 
are within practical limits (indicated below). 

It was found that the blade torsional response can be an 
essential, beneficial element of the jet flap control system. 
That is a jet-flap "torsionally controlled" rotor may be much 
more efficient than a "pure" (torsionally stiff) jet flap 
rotor where the aerodynamic control is primarily from the jet 
angle rather than the blade angle. 

Because the "mechanics" of implementing the jet-flap control 
systems were not considered, the conclusions of this study are 
independent of such practical considerations. 

The specific conclusions of this study are summarized as follows 

1) The jet-flap control system can suppress 

all the transmitted blade root shears to zero 

2) Only one independent jet-flap control mode is required 
to suppress all the transmitted shears. 

3) The jet deflection angles and additional power required 
are within practical limits. 

For example, the total jet deflection angle required 
(exclusive of the re-trim requirement) to suppress all 
the transmitted vertical and horizontal blade root 
shears never exceeded approximately 30°. Similarily 
the additional power required for the two-blade rotor 
was 117 HP, exclusive of that required to re-trim 
(301 HP). 

4. It was possible to suppress all the transmitted shears 
with the C* as low as 0.005 and the jet deflection 
angles were “still practical (at/* = 0.20, less than 
30° including requirement for re-trim). 
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5) The jet-on trim requirement adopted for this study may 
be resulting in an unnecessary power penalty for trim. 

The power required to maintain trim was greater than 
that required to suppress all the transmitted shears. 

The trim requirement adopted was that the rotor moments 
(and thus the shaft angles and fuselage attitudes) in 
addition to the forces shall be the same with the jet-on 
as with the jet-off. Accepting the small attitude 
changes may require significantly less power. The 
penalty may be due to the fact that the flapping rotor 
(with 5? 0 hinge off-set) is an inefficient moment 
generator. Thus for a semi-rigid rotor this trim 
requirement may be a reasonable one. 

6) Generally the blade dynamic bending response increased 
at all harmonics except that nearest the mode resonant 
frequency-- there it generally decreased! However, the 
net result was that peak to peak bending stresses did 
not increase significantly. 

7) Generally the more shear suppression required the greater 
the power required. 

8) Interharmonic aerodynamic coupling due to the jet may 
be quite pronounced. 

9) Suppression of one harmonic of shear affects the 
magnitude of all remaining shears--generally the shears 
at harmonics immediately above and below the suppressed 
harmonic are most affected. 

10) For the rotor studied (1st torsion at = 5.3) 

torsion is an essential and beneficial element in the 
jet flap control system. 

The blade is behaving similar to the Kaman CTR with 
the jet acting as the control flap. For the shaft driven 
rotor where the jet is heing used primarily for control 
(i.e., low Cj^ ), the blade angle-of-attack is very much 
more effective in controlling the aerodynamic forces than 
is the jet deflection angle. 

11) The power required to suppress shears and trim was 
significantly reduced by reducing the span of the jet 

by two-thirds and maintaining the same blowing coefficient. 
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7.0 RECOMMENDATIONS 


The objective of this study was to investigate the theoretical 
potential of a jet-flap control system for reducing the vertical 
and horizontal transmitted shears. It was shown that the 
potential appears good based on the computer simulation employed 
and the specific rotor configurations and problem parameters 
studied. Obviously many questions arise as to the effects on 
the results of this study of the modeling employed to represent 
the aerodynamics and structure as well as to the limiting 
assumptions made. Other questions arise relative to the range 
of specific parameters looked at; 

1) were they adequate 

2) were they representative 

3) should others be considered, etc. 

These questions obviously present valid reasons for further 
studies. 

Of all the many possible recommendations which could be made as 
a result of this study it is felt that the following are most 
important: 

1) In view of the behavior of the jet flap rotor as a 
torsionally controlled rotor, TCR, further studies to 
evaluate this possible mode of operation are 
recommended wherein the influence of the following 
parameters is investigated: 

a - Blade torsional stiffness 
b - Spanwise extent of the jet 

2) Because one of the ground rules of this study required 
that the "jet-on" rotor be trimmed( both forces and 
moments) to the corresponding values of trim jet-off, 
very large control angles and hence power were 
required just to maintain trim. This may be unduly 
penalizing the jet flap control rotor and hence the 
effect of removing the moment trim constraint should 
be investigated. 

3) Because only one calculation was made at Cj Ts( =0.005 
and the jet deflection angles required to suppress 
shear were not excessive, it is recommended that 
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further calculations be made to determine the practical 
lower bound on the jet momentum coefficient, C-Jto , and 
the jet deflection angle, H , as a function of advance 
ratio. 

4) Because no results on shear suppression were obtained 
above = 0.20, it is recommended that further effort 
be made to obtain results at a higher advance ratio. 

5) It is recommended that this investigation of the 
potential of the jet-flap for shear suppression be 
extended to include the rigid rotor configuration. 

This would provide information relative to the question 
of the possible benefits of this type of rotor con- 
figuration (with its greater control power) over the 
articulated rotor of the present study. 

6) Because the calculated power required to suppress 
shears is influenced by the value of the thrust recovery 
factor, Tr , it is recommended that calculations be made 
at additional values of Tjj • A value of 0.5 was used 
for this study, thus calculations with Tr - 0.0 and 1.0 
would establish bounds on its influence. 

7) This study has shown that the rotor blade responses 
play a very important role in suppression of shears. 

It was also noted that torsion was driven directly by 
large jet generated moments. Furthermore the large 
inherent azimuthal variations of the jet momentum 
coefficient (due only to the dynamic pressure varia- 
tions) significantly alter the jet-off aerodynamic 
mass, stiffness and damping of the blade while 

at the same time altering the motion independent 
forcing function. These observations lead to 
the speculation , that aeroelastic stability 
problems associated with such rotor systems may 
be significant. Thus investigation of such 
systems from a stability point of view is 
recommended . 
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NOTES : 0 Only vertical shear suppressed at 2nd and 5th harmonic 
respectively 

© Non-cancelling shears suppressed at 2nd through 6th 
harmonics only 

0 Short jet, only vertical shear suppressed at 5th harmonic 

0 No shears suppressed, only jet off rotor trim-forces 
and moments-required 

© Increased torsional stiffness, jet off and jet on 

(g> Convergence not obtained 

Case No. 1. 20. 03. 03 

""E-— run number within the set 

1— - des ignates tip jet momentum, Cy_ 

*o 

‘—designates advance ratio , JJ 

— designates rotor configuration 

(1) 2 bladed, 44,482 N (10,000 lb.) rotor 

(2) 4 bladed, 88,964 N (20,000 lb.) rotor 
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TABLE 4a. 

JET ON VERTICAL BLADE ROOT SHEARS - AMPLITUDE AND PHASE ( N /DEGREES) 




































































































































































































































TABLE 4/t. 

JET ON VERTICAL BLADE ROOT SHEARS - AMPLITUDE AND PHASE ( N /DEGREES) 
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JET-ON HORIZONTAL BLADE ROOT SHEARS - AMPLITUDE AND PHASE ( N /DEGREES) 
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TABLE 

HARMONICS OF TIP JET ANGLE REQUIRED TO SUPPRESS SHEARS - AMPLITUDE AND PHASE 
(DEGREES/DEGREES) 
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HARMONICS OF TIP JET ANGLE REQUIRED TO SUPPRESS SHEARS - AMPLITUDE AND PHASE 
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TOTAL POWER REQUIRED TO SUPPRESS SPECIAL SHEAR CASES 
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DESCRIPTION 
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JET ON - 0560 OHIY TO MAINTAIN JET OP* TUI* 

1 t 33 

1.20. OS. to 

JET ON - TRIM AMO SP VIST. SMEAR SUP PtESSCO . 

t Zoo 

1.20.0 VOS 

JET ON- TRIM AMO 2P VERT. SMEAR. SorpRE»€0 . 

\^\<0 

1 20. 03. 0(* 

JET ON- TRIM AMO 2Pt. <*P TRANSMITTED SMEARS 

SUPPRESS 8 0. 

1 230 

1.20.03.0% 

JET OK -TRIM AMO AIL TMKJrPTTEO SMEARS 

SUPPRESSED. 

t*3 to 

\ .20.03.S4 

> 

JET OH- SHORT JIT- TRW AMD 5? VERT. SMEAR 

SUPPRESSED . 

8 92 

2.20.00.02 

JET OPF - A BLADES ROTOR 

2 2 02, 

2.20.03.02 

JET OK - TRIM AMO All transmitted SMEARS 

-4 BLAded Rotor suppressed. 
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Figure 1. GEOMETRIC PROPERTIES OF ROTOR BLADES 


BLADE WEIGHT — 



Figure 2. BLADE WEIGHT DISTRIBUTION 
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Figure 3. BLADE PITCH INERTIA DISTRIBUTION 


EDGEWISE STIFFNESS -El/gi 



Figure 4. BLADE BENDING STIFFNESS DISTRIBUTION 
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BEAMWISE STIFFNESS -El/5x r ,f 


1-0 
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DIMENSIONLESS DISPLACEMENT — ^ /y 
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Figure 8. FIRST & SECONO TORSION - MODE SHAPE 
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FREQUENCY 




AT 


HARMONIC 



Jl = 5 

HZ 

(U/Sl) 

RIGID FLAPPING 

5.18 

HZ 

1.037 

1st FLAPWISE 

BENDING 

14.7 

HZ 

2.94 

2nd FLAPWISE 

BENDING 

26.4 

HZ 

5.28 

1st EDGEWISE 

BENDING 

5,92 

HZ 

1.18 

2nd EDGEWISE 

BENDING 

37.3 

HZ 

7.46 

1st TORSION 


26.5 

HZ 

5.30 

2nd TORSION 


69.3 

HZ 

13.86 

STIFF BLADE: 

1st TORSION 

84.1 

HZ 

16.7 


2nd TORSION 

219. 

HZ 

48.3 


Figure 9. BLADE FREQUENCIES 
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UNCOUPLED ROTATING NATURAL FREQUENCIES (cPM) 


1st TORS I ON=1 591 CPM atJl=300RpM 
2nd TORS I 0N=4l 56 CPM 



ROTATIONAL SPEED (SPM) 

Figure 10. ROTATING BLADE FREQUENCY PLOT 
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Figure 11. FORCES & MOMENTS ACTING ON THE HELICOPTER 
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ROTOR SHAFT ANGLE, -DEG 
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JET PRESSURE RATIO, VELOCITY AND COMPRESSOR POWER PLOTS 



0.20 



AZUMUTHAL STATIONS 



IX: 0.2.0 



AZIMUTHAL VARIATION OF JET MOMENTUM COEFFICIENT AT SEVERAL RADIAL STATIONS 




Figure L6. PROGRAM LOGIC FLOW 


83 





ROTOR ROLLING MOMENT - (FT -LB) ROTOR PITCHING MOMENT - (FT -LB) 
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Figure 17. ROTOR TRIM MOMENTS JET-ON & OFF vs ADVANCE RATIO 
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Figure 19. CONTROL ANGLE VALUES vs ADVANCE RATIO 
- JET-OFF 
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Figure 20. AZIMUTHAL VARIATION OF REQUIRED JET ANGLE TO SUPPRESS ALL TRANSMITTED SHEARS TO 
ZERO FOR )*. =0.08, =0.01, 0.03 



o 



ea 


AZIMUTHAL POSITION 

Figure 21. AZIMUTHAL VARIATION OF REQUIRED JET ANGLE TO SUPPRESS ALL TRANSMITTED SHEARS TO 
ZERO FOR JJl =0.20, C J T , =0.005, 0.01. 0.02. 


TRIM ©HLY 

AU TRAMS MlTTt t) SHEARS 



89 



90 


AZIMUTHAL POSITION 

Figure 23. AZIMUTHAL VARIATION OF REQUIRED JET ANGLE TO SUPPRESS ALL TRANSMITTED SHEARS TO 
ZERO FOR STIFF BLADE AT p. =0 . 20, Cj r =0.03. 





Figure 24. AZIMUTHAL VARIATION OF REQU 
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Figure 27. AMPLITUDE OF TRIM & TRANSMITTED FLAPWISE SHEARS 
vs ADVANCE RATIO, JET OFF. 
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Figure 28. AMPL 
vs 

ITUDE OF TRIM & TRANSMITTED 
ADVANCE RATIO, JET-OFF. 

CHORDWISE SHEARS 
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HARMONICS OF JET ANGLE AMPLITUDE - DEGREES 
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JET- OFF/JET- ON BLADE RESPONSES ATU = 0.20 
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Figure 31. TOTAL POWER REQUIRED TO TRIM & SUPPRESS TRANSMITTED 
SHEARS vs JET TIP MOMENTUM COEFICIENT AT = 0.20. 
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Figure 32. 


TOTAL POWER REQUIRED vs ADVANCE RATIO AT SEVERAL 
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APPENDIX I 


I. DETAILS OF THE AERODYNAMIC MODEL 


The aerodynamic model of Reference 21 was modified to incorporate 
the jet flap aerodynamics for the purposes of this study. The 
modifications were made to include the effects of the jet on the 
local blade section LIFT, DRAG and MOMENTS of a rotor blade in 
forward flight. 

The jet effects were included in the unsteady aerodynamic model of 
Reference 21 in a quasi-steady fashion. No jet dependent apparent 
mass terms were included. Similarly no explicit aerodynamic terms 
dependent upon the rate of change of the jet angle (XT ) were 
included. However, the indirect effects of the unsteady jet were 
approximately accounted for via the wake model in that changes in 
the airfoil bound vorticity are reflected as changes in the wake 
vorticity and hence modify the impressed velocity distribution on 
the airfoil. 

The following sections describe the heuristic approach which was 
employed to derive the required governing aerodynamic equations. 

Also presented are the expressions which resolve the airfoil section 
lift, drag and moment into their respective counterparts in the 
shaft plane for use in the trim equation and for determination of 
the blade root shears. Finally, the airfoil data, jet on and jet 
off, used in this study are presented. 


I. 1 Unsteady Rotor Aerodynamic Representation - No Jet 

From Reference 21, the unsteady LIFT and MOMENT are given as (on 
a given rotor blade section with no jet) 

XL = C L *. bpV,(A.+ ift.) + (3A.+A,+ 4 a») <i-i> 

'b\ = TTk?pV, ^(A.+ jA,)- 2V,H^* + 7 Ai "4 A3) J 

where Ao » A, etc. are the coefficients of the Glauert 
expansion for the bound vorticity, i.e. 

V(8,t)s zf A(t)cot(|) + ]T A»(tlsin nS^i (1 ' 3> 

M 11=1 H 


loy 


These coefficients are shown ; via the condition for no 

flow through the airfoilj to be (See Figure 1-1 for definitions) 


/ * v """ 

a. r5 v„ k = (-h +v,^ K + 1 s. K . r, 

NBA 

■ ( fc ) K - L S>KJ ^ 

HRA Ji 

= “ £ S **i H . 


(1-4) 


(1-5) 

(1-6) 


where the 


the IT. 


n K 


21 terms are the induced velocities from the wake, 
are elements of the Fourier Series describing the 


impressed velocity distribution across the chord, and 

b 


- s 


iK*) eix 


(1-7) 


the subscript K denotes the K™ collocation point in the disk; 
NRA denotes total number of collocation points in the disk. 


The relationship between the A > 5 and r*5 is obtained from 
Equations (1-3) and (1-7) with the transformation 

X = - b COS G ( I-7A) 

introduced to yield 

r K s 2 tt b K ( A. ♦ i A») K . (I_8) 

Since the theoretical circulation and lift curve slope are not 
achieved in actuality. Equation (1-8) is modified by the factor 

/2tt 

where = actual lift curve slope for a given airfoil section, 

= C L b K (Ao + -£A,) k (I- 9 ) 

Combination of Equations (1-4), (1-5) with (1-9) yields the 
simultaneous set of algebraic equations in the unknown Q which 
must be solved: 

H«A 

r K * x K + 1 <r K r] (i- io) 

j.i 


lOfc 






V,* «>s di smV* 


section a a" 


FIGURE i -I. SIGN CONVENTION USED 


IOT 



where 


C *-o<, K b k (~ h + V, ^ + 1 ^ OL^) (I-H) 

/ _ , v (1-12) 

^Kj * \ ^“ki “ z ^‘K]l 

Quasi-Steady Jet Flap 

For steady state flow, References 25 and 26 derive the LIFT 
and MOMENT coefficients for a jet flap airfoil as (respectively) 


UFT 

^obV.2 



Cl*oC -l-C. X 


' A 

Q*t 


MOMENT’ / i / A 

" ^Cj = + ^M r ”C + C M ^t 


(1-13) 

(1-14) 


where 



f s lift and moment curve slopes with 

respect to oC J Z i ‘i. 

airfoil angle of attack (2-D) (radians) 

jet angle relative to airfoil chordline (radians) 

t /2b 

maximum camber for parabolic camber (inches) 
airfoil semichord (inches) 
jet momentum coefficient 


The moment of Equation (1-14) is POSITIVE LEADING EDGE DOWN and 
is taken about the LEADING EDGE. 


The lift and moment curve slopes can be rewritten as 

Cu €l * (C^ o +C^)0C +(C«. r >S 4 (<:,.$*•*- (1-15) 

C Htfj * ( 4- (C^ X )Z + + (1-16) 

where 


C Uo( , C , etc. are the airfoil lift and moment curve 
* slopes with no jet 



C ^ * C^, etc. are the jet on airfoil lift and moment 
curve slopes as functions of Cj only 

But Equation (1-16) is for a positive, leading edge down moment, 
and must be transferred to a positive, leading edge u£ moment 
about the airfoil midchord to be consistent with RAS (See Figure 
1-1). Thus, 


C He J 


MOMENT ABOUT THE MIDCHORD, LEADING EDGE UP 
POSITIVE 




(I-l 7 ) 


Now introducing the quasi-steady approximation for the effective 
angle of attack , 

<* = - 7 ( «- l8 > 

o 

and recognizing that, quasi-statically, a pitching rate, OCo , 
is equivalent to a parabolic camber of 




t e = 


4V, 


(1-19) 


with positive, leading edge up , 

then the total effective camber is given by 

A A A 

t =t. *t e 

A 

where - t # = built in parabolic camber, 


(1-20) 


and Equations (1-15) and (1-17), upon substitution of equations 
(1-18) and (1-20) into them, become 



♦&*- roK-5, ) ♦ (c« t - iw + & 


( 1 - 21 ) 


( 1 - 22 ) 


1.3 Combination of Jet Flap Quasi-Steady and Rotor Unsteady Aerodynamics 


If a comparison is made of the quasi-steady terms of Equation 
(I-l), i.e. 




(1-23) 


* Assuming uncambered airfoils, 



with the quasi-steady approximation for the jet flap in Equation 
(1-21), appropriately modified by ( ), it is recognized 

that the non- jet dependent terms are readily identifiable (Note 
that Ci$ # « dtifbt, = 4* or applying same correction on the 
theoretical lift curve slope ). Thus if the 

quasi-steady terms of Equation (1-23) have added to them the 
remaining jet dependent terms of Equation (1*21) the resultant 
equations would represent the quasi-steady aerodynamics of a 
jet flap rotor in a "strip” sense. The remaining wake and un- 
steady terms are now included to complete the aerodynamic 
representation for the non-jet model and also allows an approx- 
imation to the jet dependent wake effects. This occurs due to 
the fact that the strength of the wake vortices will be dependent 
upon the blade bound vorticity which is proportional to the local 
load which is directly dependent upon the jet via Equation (1-21). 

Hence the resultant LIFT expression for the jet flap rotor becomes 

£ ■ (a.+ jA.) tWf i(}A.+A, + ^A,) 

+ Cl* bfV, (*•") + Ci^b^'M + CLfcbpV, (;|A,) ' 


A similar arguement for the moment yields 

V\.= nb* f V.J(A.*iA.') - ^i(A.4»A..^A»')|] 

-2 f bV,<C M ,-iCm)(A.) - 


- ?^bV, (C M $ . 


(1-25) 


Since the lift must be proportional to the total bound circula- 
tion plus the noncirculatory terms in the jet flapped airfoil 
case as in the classical case (Reference 24), the relation may 
be written 


= f.V.r + 7Tb l ? + 

which when compared to Equation (1-24) leads to 

r K «| C L „ o b(A.+iA. > ) + C t- b(A.) 

4- Ct x loV ( iX ) + b (^A J 


(1-26) 


(1-27) 


no 



Equation (1-27) may be rewritten to yield a form identical to 
Equation ( I - 1 0 ) 

K3A 

& ' Ik + Z R (I ‘ 28> 

but now where * 

' Ct *o K b K (“^ + 

+^«t«b K (-i+v,o^) + Ct tH bV lK rK (I . 29) 

. (I_30) 

Thus Equations (1-24), (1-25) and (1-28) when coupled with the 
expression for the drag at each disk collocation point completely 
define the aerodynamics of the rotor for the model assumed. 


1.4 Section Drag Expressions 


The airfoil section drag, © -as defined in Figure 1-1, was 
described as composed of two terms. The first term represented 
the jet off blade profile drag and was assumed to be a function 
of the freestream dynamic pressure, the local mach number ( M ) , 
the local effective angle of attack and the airfoil profile. It 
was assumed that this term was unaffected by the jet. The second 
term represented the direct jet momentum reaction parallel to the 
blade chord line and a term describing the jet momentum recovery. 
Thus, the drag expression was written as 


(1-31) 


The jet dependent term, , is defined as 

T r +0-Tr)cO 5(«^+'O ( l /r+JJlSin'A) (1-32) 


where 


Tr = 


thrust recovery constant 


= total 


loca 1 


instantaneous geometric angle of attack. 


1.5 Shaft Plane Forces 


Resolution of the local airfoil aerodynamic forces and moments 
into the forces and moments (locally) referenced to the shaft 



plane in the rotating frame was accomplished via the following 
expressions (See Figure 1-1) 


f; 

m 

e£smfl(i + ©co»<^£ 

(1-33) 

Fy 

ZB 

1 a? 1 hi ) 

(1-34) 

Ft 

- 

adcesotl - © smoCi. 

(1-35) 



(1-36) 

Under the 

assumption of small angles ( dl«i ), Equations 


(1-33) to 

(I- 

•36) become 


Fjc 

a 

X + © 

(1-37) 

f; 

a 


(1-38) 

F ; 

- 

X -®(*0 

(1-39) 



(1-40) 


The local induced angle oti , was written as 


oL{ - — x 

V, 


TOTfct neUTP/C A« VEL OCITY HOttMAl shaft plame 
~ 


where 


Vj sm oil * V^cosots cos^ (ji e ♦ 1 37* ^ ) 

+!*„>; -cltjCj + vr, j 


(1-41) 


p t - precone angle 

4^. - mode shape of i ch flatwise mode 
hj - tip deflection of i cb flatwise mode 
4^ - mode shape of 1 tb torsional mode 

- tip deflection of < cb torsional mode 
4 t . - mode shape of j cb control mode 
Cj - tip deflection of j th control mode 
€ - distance from elastic axis to midchord 


(IZ 



®o - distance from control axis to midchord 
- induced velocity. 

The forces and moments given by Equations (1-37) through (1-40) 
were then used to determine the generalized forces and hence the 
blade root shears. 


1.6 Airfoil Character is tics 

The lift, drag and moment expressions of Equations (1-24), (1-31) 
and (1-25) respectively were defined in terms of jet off and jet 
on airfoil characteristics, 

(Cm* } 

and 

( 

These airfoil data are required as input to the RAS program; 
theoretical or experimental values may be used. For the jet off 
coefficients, Ci^. , and a radial distribution 

may be inputed. These values are used up to the stall angle. For 
the present study the following values were used as obtained for 
the NACA 0012 airfoil from Reference 28 

=5.73 

» 2 Cl*. 

Ch*." ( 

■ c »-*. 

The corresponding experimental airfoil drag data (Reference 28) 
required a two-variable polynomial fit. The variables were 
otm f H ( M = Mach no.). A third order polynomial on <*e was 
fit for the Cp^ vs o<* curve at four discrete mach numbers (0.45, 
0.50, 0.60 and 0.80). The resulting set of coefficients were 
then cross fit as a function of M . Thus 

c 0p s ao<«)+a,tH)o(e4ajW«| +a%wotl 

For each M , a value of ofumr was selected at which the value 
of Cp r was limited. For all o< e above ^ ltw • 



This was done primarily to simplify the polynomial fits required 
to represent the variation. A plot of the computed vs Ote 

for several M is presented in Figure 1-2. 

The jet on coefficients, Cm* etc. were adapted from Hough 

(Reference 25) where plots of etc. vs Cj are presented. 

Polynomial fits on Cj were made for each aerodynamic coefficient. 
A synopsis of the polynomials used in this study along with the 
range of Cj applicability of each is presented in Table 1-1. 

Also presented are the jet off values. Presented in Figure (1-3) 
are plots of these jet on coefficients for the range of Cj en- 
countered in this study. Note that Cw t and £*$ are plotted as 
their negative values. 


1.7 Treatment of Stall and Reverse Flow 


The stall characteristics of the jet flap airfoil have been 
assumed to be substantially those of the jet off airfoil. Thus 
the stall angle for the no jet airfoil will remain the stall 
angle for the jet on airfoil. However the maximum attainable 
lift coefficient will be larger by an amount proportional to 

C T > • 


An airfoil section is defined to be stalled when 


or 


where 


°*e = , for *V, 

>0 iL , for - (reversed flow) 

V, 


! 


input stall angles for forward and reversed 
flow respectively 


(1-42) 

(1-43) 


A» uniform component of chordwise distribution of 
impressed velocity as per (1-4). 


The expressions defining the lift, drag, pitching moment, and 
total bound circulation of the airfoil section operating above 
the stall angle of attack are presented below. 

For the total bound circulation, a maximum or limiting value in 
defined based on Od H , thereby 

* tV| K £ (Cc* # + Ctoc) Z ( 1-44) 

where =<* H or 0 ^ depending on the sign of . 
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igure I -2 PLOT OF nRAG COEFFICIENT vs ANGLE OF ATTACK 


TABLE 1-1 


SPECIFIC EXPRESSIONS FOR AIRFOIL COEFFICIENTS 
AND RANGE OF APPLICABILITY 


SLOPS 

JIT 

*MOtP. 

JIT DIPPNOONT 

Cj RAH*f 

Cu* 


C^(«WI c/* 4 0*20 Cj ) 

O^CyiS-O 


o 

(2C^C,) (l+ 0.181 c/ 1 4 0.15* C,) * 

0*Cj$S.O 

Ci t 


(l«t«.)(l.*»*c v vm Cj'-o.'so^C] .amcj) 

H 

C H* 


( c k./<X°»’ c J-° 0|S ’ C J ? 

■ 

C WT 

0 

Ci*. (<UI« Cj + 0.2*J Cj* ) 

0* Cj<S.O 

Cwf 

Cw. 

Ct*. (oWS C 3 - O.THCj l 4 0.010 cj - 0. Olfcicj ) 
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LIFT CURVE SLOPE COEFFICIENT MOMENT CURVE SLOPE COEFFICIENT 



For a small range of ole from o£m to oL , it was assumed 
that T* remained constant and equal to Thk ; for oie from SL 
to Tt/z f r K , was assumed to vary linearly to zero. The 
defining equations are 


Oi » ( 1.25) 
for oi ) 

for ( 5<l*«»Ot/2 ) 


r « , ( r «-(^a) (w ‘ uS ‘ ) ] (SaiH “'* ) 


(1-45) 


(1-46) 


(1-47) 


A pictorial description of the above is presented in Figure 
(1-4). 



Figure 1-4 STALLED GAMMA CHARACTERISTICS 

Thus during the solution the magnitudes of the T which are shed 
into the wake are limited to reasonable values. This is an 
approximation to prevent large, unrealistic values of vorticity 


from being shed into the wake from sections which are above their 
stalling angle of attack. 


The lift and pitching moment are computed by (1-24) and (1-25) 
when the airfoil section is below stall. But when the section 
is stalled, the lift is approximated by the following: 


* P^»k Iw* + A % 


(1-48) 




where 


Ai k r 2^b K Vt k stn 1 o<t K 


That is, the stalled lift is computed as the sum of the stalled 
(limited) circulatory lift and a cross-flow drag force, A£ k • 


The corresponding stalled moment is computed as 


*1 


Mk 


where 


CP * 


oC Mk m ^ 

V> CN«I ~ *1 


(1-49) 


Thus the stalled pitching moment about the midchord is computed 
as the moment due to the stalled lift, 3 £mk acting at a center of 
pressure which moves toward the midchord as o<e approaches IT fa 
At oftrTT/l , the pitching moment goes to zero as might be expected 
for a flat plate normal to the direction of flow. The velocity, 
Vt , is the total local resultant velocity acting normal to the 
spanwise axis (or midchord axis) of the local blade segment and 
is given by 


t 


( t l -v v r v 


(1-50) 


The velocities V, and Vj, are the two orthogonal components of 
V T in a plane normal to the local spanwise axis. Velocity V| 
lies on the intersection of this plane with a plane normal to 
the shaft. 

In reverse flow regions, jet off, the same approach as outlined 
above and depicted in Figure 1-4 (see reverse flow region) is 
employed with appropriate sign changes to define the blade 
sectional forces and moments. 


With the jet on in reverse flow it was assumed that the effect 
of the jet would be to cause flow separation and hence the airfoil 



was assumed to be stalled for all o U • Thus the corresponding 
lift and moment were computed as 

C, K smT* (1-51) 

(jet on, reverse flow) 

H - C, R smT K d-52) 

The expression for the drag given by Equation (1-31) includes 
the influence of stall through the Cop characteristics used 
(Figure 1-2). The direct jet momentum contribution to the drag 
is assumed to be unchanged by stall. 



APPENDIX II 


II. 0 ROTOR TRIM EQUATIONS 

The rotor trim equations incorporated into the formulation 
are for the zeroth harmonic rotor forces and moments as 
defined in the non-rotating shaft oriented coordinate system. 

The equations for these trim parameters are written in terms 
of all six components of the blade root forces and moments. 

Thus these equations are explicit functions of the blade 
response, the aerodynamic loads, and the rotor control 
. variables (i.e., all the other dependent variables). They 
are based on the development by Chang, presented in Reference 
29, and were generalized to accommodate most any combination 
of blade root constraint. They contain non-linear terms which 
depend on steady-state values of parameters such as flapping 
and lagging displacements and, for a blade without a lag hinge, 
the steady dragwise root moment. Because these steady-state 
parameters are constants which can readily be estimated, the 
non-linear terms are accommodated by inputing and using the 
estimated values. 

The different blade root constraints give rise to different 
terms, thus the form of the equations change with root 
constraint. For example, a blade without flapping hinges 
has a large hub moment contribution from the lift loading, llr) 
given by, ? 

= jrloridr 

Whereas a blade with a flap hinge has, instead, a term given 
by K 

Mj =■ r f 

0 

where ** is the flap hinge offset from rotor shaft. Many 
similar situations exist for the rotor forces and moments due 
to blade inertia, centrifugal, and aerodynamic loadings. 

In the computer program, the required change in the form of 
the trim equations with blade root articulation is implemented 
by including in the equations the terms for all the different 
root conditions. Each of the terms have additional coefficients 
which are either zero or one depending on the articulation. 

Thus the appropriate terms (or parts thereof) are kept and/or 
deleted in various combinations as required. This implementa- 
tion requires only two inputs, F and L.. These inputs 
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represent the presence or absence of flapping and lagging 
articulation respectively by a value of 1 or 0. 

There are three rotor trim forces and three moments--all are 
included in the analysis. Of the forces, one component (that 
normal to the shaft-plane) is already represented by the sum of 
the zeroth harmonics of the blade root shear, so that it was 
only necessary to add the equations for the remaining two inplane- 
components of rotor trim force. Relative to the moments, because 
only steady trimmed flight is considered, one of the moment 
equations (for yawing moment) is implicitly satisfied (eliminated) 
by always constraining the rotor side force to be equal to the 
value of the rotor torque divided by the tail rotor moment arm. 
Thus only the equations for the remaining two rotor trim moments 
(pitch and roll) need be added to the analyses, Therefore of the 
six force/moment equations describing the rotor trim, only four 
had to be added to the analysis. 

Each of these four trim parameters is primarily a function of 
the first harmonic blade responses and airloads, thus their 
equations are included in the first harmonic equation-set. The 
equations presented on the following pages are in the form as 
they would appear in the set of simultaneous equations at the 
first harmonic. In these trim equations FNY\ and L*A\ defined as 

F(vM * (F 

are computed internally and used in addition to P and L 
to include and delete terms according to the articulation. 

The remaining symbols represent physical parameters of the 
rotor and blades, and the harmonic components of the blade 
responses, and airloadings as defined at the end of this 
Appendix. The dependent variables of the equations are 
underlined. 
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PITCHING MOMENT TRIM EQUATION : 
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ROLLING MOMENT TRIM EQUATION : 
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LONGITUDINAL FORCE TRIM EQUATION ; 



LATERAL FORCE TRIM EQUATION: 
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Nomenclature for Trim Equations 


- Built in precone 

^ } Q - Offset from rotor shaft of flap, laig, 

and equivalent precone hinges 

R 



o 


spanwise mass distribution 

flapwise mode shapes 

edgewise mode shapes 
torsion and pitch mode shapes 

torsion or pitch mass moment of 
inertia distribution 


where WV 



[ 37 + &] ' 



The following are harmonic components of the dependent 
variables (blade response) 

* 1st harmonic flapping, cosine and sine components 


(ft. hj)) 
tft.Hi) , 
(ft. 9;) 
(ft.O, 


1st harmonic cosine components flapwise, edgewise, 
torsion, and rigid body pitch 


( &.hi)) 

(6.HjU 

( 6.6j) 

(6,0 


1st harmonic sine components of the above 


c' c' c ' 


zeroth, and the 1st harmonic cosine and sine components 

of R 

r F 

where is the aerodynamic drag distribution 


The following are generalized aerodynamic forces: 



1st harmonic cosine 

where the Ftr'l 
for ho , Ho , C, . 


components of 

is t(r\ ,£(r\WKO 


respectively 


1st harmonic sine components of the above 
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APPENDIX III 


III.O JET COMPRESSOR POWER REQUIRED 

The total power required for the jet flap rotors studied 
herein may be written as 

P T ■ Pr. + P c (III-l) 

where 


= shaft power required by the Rotor^HP, 

R: = compressor power required to supply the jet, HP. 

The rotor shaft power required, P R , is computed directly as 
an integral part of RAS program as 


P R , Nil 

Ztt (550) 


R 


where 


O o 


Arif 


(ni-2) 


N =■ number of blades 

(see Eq. 1-37 of Appendix I). 


To determine the compressor power, P c , a simplified analysis 
for a one -dimensional, compressible channel flow was employed. 
It can be shown that for isentropic flow 


K - 


1*V 


a-R*T«M l| (y- /tw V .|J] 


(III-3) 


where 


5S0U-I) nt 

total Rotor mass flow rate (i.e., for 
all blades 

"P, - total pressure at compressor outlet 

"P - atmospheric pressure 
ATM 

W atmospheric temperature 
R - gas constant 

- ratio of specific heats, (1.4 for air) 



“Y\ c - compressor efficiency (assumed at 0.8) 

2 

- gravitational acceleration, (32.2 ft/sec )# 

There are two unknown parameters in Equation (I1I-3)* they are 
the total Rotor mass flow rate for all blades, Mj 7 and the 
pressure ratio (^. /-p*** ). 

• 

The mass flow rate, , may be determined as a function of 

the pressure ratio by considering the mass flow rate per unit 
span from a converging tube or orifice for an adiabatic 
frictionless flow (see Figure ITtl). 



Figure III-l ASSUMED JET DUCT CONFIGURATION 


The mass flow per unit span is thus given by 

where hj = height of the jet slot, ft 
t J ft/ see. 11 

•• . . 

= “P* - (duct losses). 

The total mass flow rate required of the compressor for N 
blades is then given by 


r R 

Mj t » Nj ^ <jpj ^ 


(III-5) 


12.9 



where Tj = inboard radial extent of the jet. 


If the radial distribution of the jet slot height, 
given by 





TlF 


*2 


is 


(III-6) 


then the expression for the total rotor mass flow rate as a 
function of the pressure ratio is obtained from Equation III-5 
upon substitution of (III-4) and (III-6), 



(III-7) 


An expression for the pressure ratio, ( Ip. /'f ATn ), can be 
obtained by considering the jet velocity, Vj , at the exit 
of the nozzle of Figure III-l, 



( III-8) 


If the jet velocity is required to be directly proportional to 
the local freestream velocity, it can be written as 


Vj » fe. (Ril)( r +jUL3 ln ^) (111*9) 

where k» is the proportionality constant. 


Then using Equation (III-9) in Equation (III-8) and solving for 
( IP. /fvv ) an approximate expression for the pressure ratio 
is obtained as 




+ 2k* (R-fi} C.r^u-suh'/' 

i (m-io) 


-f* k. (Ril) C* ( |xsm i> ) 


i h 


where C. - 


i - 1 


(III-10A) 


Thus from Equation (III-10), ( 1 1 1— 7 ) and ( III- 3 ) , the compressor 
power may be determined as a function of ambient air condition, 
advance ratio, the slot height distribution and duct pressure 
losses. 



For the cases analyzed, the following siroplifing assumptions 
were made: 


(i) the azimuthal and spanwise distributions of 
the jet velocity are constant and equal to 
the maximum rotor tip velocity, i.e. 

» k. ( &XL)(l+^) 

(ii) centrifugal pumping action was just 
sufficient to balance duct pressure 
losses and hence 

Ip. • T* 

(iii) the atmospheric condition to which the jet 
exhausted were those of the ambient air- 
stream at infinity. 


(III-ll) 

(III-11A) 


The assumption of Equation (III-ll) yields for Equation (III-10) 




(RXlf C + 2 km C -° JA- 



(WO Co^L 

\ 


. (III-10B) 


The jet slot height at 
of C 3% as 



the tip may be written as a function 


C Jt I 1 

k.( •+/»■) i 



(III-12) 


which when substituted in the Equation 

W T r C T N(bfRil') 

J T J T- 


(III-7) yields 



( 1 1 1 - 1 3 ) 


Presented in Figure 13 are plots of Vj , 1 »\ 3t , and T, 

vs ju. for various . The following are the values of the 

parameters employed in this study. 


N = 
R-ft. ■ 
R = 

i = 

b = 


2 


- 1.2 

760.3 ft/sec. 


= 12.1 ft 

24.2 ft 

f 

= 0.002255 slugs/ft 

1.4 

R 

= 53.3 ft/ a Rankins 

0.755 ft 
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APPENDIX IV 


IV. 0 DETAILS OF THE PROBLEM FORMULATION 


Fundamentally the overall problem is formulated in terms of 
more (two) dependent variables than equations. Two constraint 
equations are then added to this equation-set so as to make the 
number of equations equal to the number of variables (effective- 
ly reducing the number of variables by two). 

The nature of the problem is such that the solutions are 
periodic on the rotor speed; thus expanding the solution in 
a Fourier series (Section 4.2) results in doubling the number 
of equations and dependent variables. That is each basic 
dependent variable of the problem is represented by its 
amplitude and phase (sine and cosine components). With the 
exception of the rotor trim parameters, all dependent variables 
are defined in the rotating shaft oriented coordinate system. 

The rotor trim parameters are defined in the non-rotating shaft 
oriented coordinate system. 

The formulation is implemented such that the amplitude and 
phase of any two of the dependent variables may be constrained 
individually at each harmonic. The following paragraphs des- 
cribe this overall problem formulation in more detail. 

Classically the set of equations describing a system is written 
such that all the terms which are functions of the dependent 
variables are grouped on the left-hand side of the equations. 

The remaining terms (not functions of the dependent variables) 
are then grouped on the right-hand side of the equation and 
are commonly referred to as the forcing functions of the system. 
However it is not possible to formulate the aeroelastic response 
problem in this manner with sufficient accuracy, as explained 
in the following. 

The rotor aeroelastic response problem can be considered as 
being composed of two separate but interacting aspects -- they 
are a structural dynamics representation and a representation 
of the system aerodynamics. In reality these two aspects are 
intimately interrelated. Because it is not possible to accurately 
express the aerodynamic forces experienced by the blades as 
explicit functions of the blade response (including the effects 
of the wake, etc.), the aerodynamic forces are computed by what 
is more properly described as a numerical simulation of the 
aerodynamic aspect of the problem. That is, given the rotor 
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operating conditions and the blade response, the simulation 
computes the resulting aerodynamic forces as functions of time 
such that they are consistent with the blade response, the wake 
of the individual blades, and the non-uniform inflow (as dis- 
cussed in Section 4.1). These resulting aerodynamic forces 
are, however, the sum of the response dependent airloads and 
the response independent airloads (i.e., the forcing functions 
for the dynamic system). Thus the aerodynamic forcing functions 
and the response dependent aerodynamic forces are not available 
separately, nor is a description of the latter available as 
explicit functions of the dependent variables for writing the 
equations. That is, with a numerical simulation of the aero- 
dynamic aspect of the problem, it is not possible to write the 
equations of the system in the classical form. The procedure 
for handling this apparent difficulty is part of the solution 
procedure and is described in Section 4.4. 

The vertical and inplane blade root shears are functions of the 
blade dynamic response as well as the aerodynamic loadings and 
thus were treated as two additional basic dependent variables. 
That is, the equations for these two components of blade root 
shear are added to the seven equations of blade motion at each 
harmonic (except the first). 

Because the longitudinal and lateral inplane rotor trim forces 
(i.e., the zeroth harmonic components in the non-rotating 
reference system) are determined by the first harmonic blade 
root shears, these shears cannot be constrained independent of 
the rotor trim forces. Thus, at the first harmonic, the rotor 
trim equations (forces and moments) were added to the seven 
blade equations of motion, instead of the blade root shear 
equations, because these trim parameters depend on the blade 
dynamic response and aerodynamic loadings, at the first 
harmonic. At the zeroth harmonic the vertical blade root 
shears are the rotor force normal to the shaft plane so no 
modification was required to include this trim parameter. 

Included in each of the resulting nine equations at each harmonic 
are terms which explicitly depend on each of two control para- 
meters. They are the rigid blade pitch angle (collective and 
cyclic) and the jet flap deflection angle. Associated with the 
rigid blade pitch variable are mass, elastic, centrifugal, 
gyroscopic, and aerodynamic coupling terms in each of the nine 
equations. The jet flap variable has no mass/elastic character- 
istics and therefore has only aerodynamic coupling terms in each 
of the nine equations. These two control variables are not 
considered as "degrees-of-freedom", i.e., there are no equations 
of motion representing them. However, they are treated as 



dependent variables in the equation set. Thus there are eleven 
basic dependent variables but only nine equations. 

Two "constraint" equations are, therefore, added to the equation 
set (resulting in a set of eleven equations) which in effect 
reduces the number of basic dependent variables to nine so that 
the number of equations and variables are equal. As explained 
below, these two constraint equations are used to select which 
and to what value two of the eleven basic dependent variables 
are to be constrained. 

All of the coefficients in each of the constraint equations are 
zero except for a "one" at the "column-position" of the variable 
to be constrained. The right hand side of each constraint 
equation is set equal to the desired constraint value for the 
constrained variable. Thus, for example, if eauations seven and 
eight of the eleven equations (presented here in matrix 
notation) . 


(4. 3.-1) i *••»« 

are the constraint equations and it is desired to constrain, 
say, x* and it, respectively to the values A and B, then the 
"one" would be put in column positions 4 and 9 of Equations 
7 and 8 respectively and set equal to A and ;j t » B. 

Thus Equations 7 and 8 of (4.3-1) would be 

O'X ) +0'K t +0>)! 1 +l > ^+O'Hs + ••• + O-Xh * y, - Pt 
0-X, * O.-Xv * • • 0‘)6„ * - $ 

which reduce to 

X A » ft 
*¥■«* =■ 8 . 

The variables selected to be constrained (i.e., the positions of 
the "one" in each constraint equation) and the constraint values are 
controlled by program inputs. 

Because the problem is periodic and the number of equations is 
twice the number of basic dependent variables (as explained at 
the beginning), there are actually 22 simultaneous equations at 
each harmonic for the 11 dependent variables. This quite general 
formulation of the aeroelastic response problem thus allows the 
amplitude and phase of any two of the eleven dependent variables 



at each harmonic to be constrained to any desired value and the 
resulting solution obtained. It can be used to perform the 
direct rotor aeroelastic solution by constraining the pitch 
control at IP, i.e., given the operating condition and control 
setting, compute the resulting unsteady airloads, blade 
responses, rotor loads and moments, and performance. It can 
also be used to determine the pitch control required to attain 
trim by constraining the rotor trim forces and moments at IP. . 
Or it can be used to perform either one of the above solutions 
with constraints impressed at the higher harmonics on any of 
the variables. Thus the method has great versatility and 
power . 



